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ABSTRACT 

We present mid-infrared spectra and photometry from the Infrared Spectrograph on the Spitzer 
Space Telescope for 51 OH megamasers (OHMs), along with 15 galaxies confirmed to have no 
megamaser emission above Lqh = 10^'^ Lq. The majority of galaxies display moderate-to-deep 
9.7 /im amorphous silicate absorption, with OHM galaxies showing stronger average absorption and 
steeper 20-30 /im continuum emission than non-masing galaxies. Emission from multiple polycyclic 
aromatic hydrocarbons (PAHs), especially at 6.2, 7.7, and 11.3 /im, is detected in almost all systems. 
Fine-structure atomic emission (including [Ne II], [Ne III], [S III], and [S IV]) and multiple H2 rota- 
tional transitions are observed in more than 90% of the sample. A subset of galaxies show emission 
from rarer atomic lines, such as [Ne V], [O IV], and [Fe II]. 50% of the OHMs show absorption from 
water ice and hydrogenated amorphous carbon grains (HACs), while absorption features from CO2, 
HCN, C2H2, and crystalline silicates are also seen in several OHMs. Column densities of OH derived 
from 34.6 /im OH absorption are similar to those derived from 1667 MHz OH absorption in non-masing 
galaxies, indicating that the abundance of masing molecules is similar for both samples. This data 
paper presents full mid-infrared spectra for each galaxy, along with measurements of line fluxes and 
equivalent widths, absorption feature depths, and spectral indices. 
Subject headings: masers - galaxies: interactions - galaxies: nuclei - infrared: galaxies 



1. INTRODUCTION 

OH megamasers (OHMs) are 18-cm masers with in- 
tegrated line luminosities on the order of 10^ — lO'' L0. 
They are an extremely rare phenomenon in the local uni- 
verse, with roughly one hundred current ly known out to 
a redshift of z = 0.265 ()Baan et al.llT99^ . Ah OHMs, in- 
cluding the more powerful "gigamasers" {Lqh > 10'* L©) 
are associated with starburst nuclei in merging galax- 
ies. OHMs have been identified in many different types 
of nuclear environments as classified by optical spec- 
tra, but the merging galaxies are without exception (ul- 
tra)luminous infrared galaxies ([U]LIRGs). Since OHMs 
are signposts of gas-rich merging galaxies, their presence 
can also indicate the existence of associated phenomena 
including massive black ho le mergers and highly obscured 
circumnuclear starbursts (jParlingl |2007| ). OHMs are a 
powerful tool in this respect due in large part to their 
ability to be seen at cosmic distances. In order to employ 
OHMs as tracers, however, the assumption must be made 
that the OH line properties remain constant as a function 
of cosmic time and host environment. An explanation of 
the physical mechanisms and conditions responsible for 
distinguishing OHMs from non-masing ULIRGs is thus 
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vital for understanding both the megamaser phenomenon 
and the associated merger characteristics. 

Spectroscopic studies of the mid-IR emission in the 
host galaxies offer multiple diagnostics which can pro- 
vide clues to the nature of the maser pumping mechanism 
and the associated OH emission. We used the Infrared 
Spectrograph (IRS) aboard the Spitzer Space Telescope 
(jWerner et al.ll2004( ) to study merging ULIRGs. Since 
the dusty nuclear regions are typically obscured at opti- 
cal wavelengths, mid-IR observations can yield valuable 
information specific to the locations in which the OHMs 
are generated. These include measurements of the dust 
temperature and optical depth (from broadband photom- 
etry and absorption features), the excitation and tem- 
perature of the gas (molecular and fine-structure atomic 
lines) , and high-ionization lines that can signal the pres- 
ence of an active galactic nucleus ( AGN) , a possible heat- 
ing source for the dust. 

This data paper presents full low- and high-resolution 
IRS spectra along with measured mid-IR properties for 
51 OHMs and 15 non-masing ULIRGs. An accompany- 
ing paper (Willett et al. 2011; Paper II) presents the full 
analysis, statistical comparisons of the masing and non- 
masing galaxies, and tests the viability of current OHM 
pumping models based on the IRS data. 

2. THE SAMPLE 

The OHM host galaxies selected for IRS observations 
were primarily drawn from the Arecibo OHM survey 
(iDa rling fc Giovanellill2002al lbl). We selected well-studied 
OHMs with unambiguous maser detections and large 
amounts of ancillary data (including OH line and radio 
continuum maps, near-IR imaging, and optical imaging 
and spectroscopy) to maximize scientific return on the 
sample. A lower threshold of Lqh > 10*'^ Lq also elimi- 
nated extragalactic "kilomasers" from the sample, which 



2 



WILLETT ET AL. 



are likely po wered by different radiative processes than 
megamasers (IHenkel fc Wilsonl [1990!) . 

In order to be detected in reasonable integration times 
using the IRS, we required that all potential targets 
have S{60 ^m) > 0.8 Jy as measured by the IRAS 
satellite. After removing objects already observed by 
the IRS (largely th rough the GTO ULIRG program; eg 
lArmus et al.ll200"7l ). we observed 24 galaxies in the red- 



shift range 0.1 < z < 0.2. 

We supplemented these galaxies with additional spec- 
tra of OHM hosts publicly available through the Spitzer 
archive. To ensure uniformity of the data, we selected 
only galaxies from the archive that had full coverage with 
both the IRS low- and high-resolution modules. As of 
March 2008, the publicly available data from the archive 
yielded an additional 27 OHM galaxies (Table [T|). 

In order to provide a baseline for analysis of the OHMs, 
we also identified a control sample of ULIRGs that 
showed no megamaser emission above a firm limit. To 
identify these galaxies, we drew on non-detection s from 
OH su r vevs bv | Baan et ahl (. 1992): Stavclcv-Smit h^t al.l 
Jl992); Darling & Giovanellii (|200a l200lL l2002al ) and 
|Kent et al. (2002). The upper limit for OH emission is 
conservatively derived from the rms noise in the spec- 
trum at 1667 MHz, assuming a boxcar line profile with 
a linewidth Av — 150 km/s and a 1.5(J detection: 



T max 
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For this control sample, we set an upper limit of 
< 10^'^ Lq; this limit compromises between en- 
suring that all but the faintest megamaser emission is 
excluded and yielding a reasonable number of objects in 
the control sample for statistical analysis. All 51 OHMs 
in our IRS sample have Lqh above this limit. 

In addition to selecting galaxies based on OH non- 
detection, we imposed two additional criteria to en- 
sure that the control sample was as similar as possi- 
ble to the OHM hosts. Firs tly, we set a lower limit 
on the far-infrared luminosity (iSan ders fc Mirabeilll996[) 
of the non-masing galaxies as measured by their IRAS 
fluxes. OHMs occur exclusively in IR-bright galax- 
ies, due to the fact that the maser is pumped primar- 
ily by rotational transitions of a few hundred K abov e 
the ground state jB aan et al.| [T98l IHenkel et al.lll987[ ). 
iDarling fc Giovane lli (2002a) show that the relationship 
between the OH and infrared luminosities is a power law 
with Lqh Lpjjj^- Since no OHM observed with the 
IRS has LpiB. < IO^^Lq, we established this as the lower 
limit for inclusion in the non-masing sample. 

Secondly, a cutoff in redshift space is applied to sam- 
ple a sufficiently large volume (V ^ 1 Gpc^) in or- 
der to avoid systematic effects such as the Malmquist 
bias. The available data in the archive contained many 
more galaxies at lower redshifts {z < 0.05) than those 
further away. To avoid over-weighting the control sam- 
ple towards galaxies at low redshifts, we sorted galax- 
ies that met the Lq^^ and Ljj^ criteria into bins of 
A z — 0.02. For bins where the number of non-masing 
galaxies exceeded those of OHMs, we randomly removed 
objects from the non-masing bins until the numbers were 
equal. This reduced the control sample to 1 galaxy at 
< z < 0.02 and 4 galaxies at 0.02 < z < 0.04; for aU 
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Fig. 1. — Distribution of integrated OH luminosity for the IRS 
samples as a function of redshift. For non-masing galaxies, the 
OH luminosity is an upper limit calculated from Equation [T] The 
dashed line at Lqh = 10^'^ Lq is the upper limit on possible OH 
emission for the non-masing control sample. 

other redshift bins, no such adjustments were necessary. 
As of March 2008, there existed 15 suitable candidates 
in the Spitzer archive qualifying for the non-masing con- 
trol sample (Table [1]). Figure [1] shows the distribution of 
OH luminosity for all objects as a function of redshift; 
for the non-masing galaxies, we display upper limits as 
computed in Equation [TJ 

Throughout this paper, we assume the WMAP5 cos- 
mology with Hp = 7 0.5, Qm = 0.274, and 17a = 0.726 
(jHinshaw et al.ll2009l ). 

3. OBSERVATIONS 

We observed 24 OHMs with the IRS from August 
2006 through December 2007. The IRS contains four 
modules in two different spectr oscopic resolutions : low- 
(LR) and high-resolution (HR) ()Houck et al.|[200l . The 
short-high (SH) and long-high (LH) modules operate at 
a resolution of i? ~ 600; the short-low (SLl fc SL2) and 
long-low (LLl fc LL2) modules operate at resolutions 
ranging from i? ~ 56 — 127, depending on the observing 
wavelength. Two orders, SL3 (7.3 — 8.7 fim) and LL3 
(19.4 — 21.7 fJ-m), cover the overlapping range between 
the first and second orders in both SL and LL. We used 
these only as checks for the absolute fiux calibration be- 
tween the different orders. 

All targets used the Staring Mode Astronomical Ob- 
serving Template (AOT) with the galaxies placed at two 
nod positions approximately one-third and two-thirds 
the length of the slit. We observed targets in all six 
modules as well as an equal-time, off-target exposure in 
the LH module to be used for sky subtraction. 

Although the majority of our sources had no previous 
IRAS detections in either the 12 or 25 /xm bands, we ex- 
trapolated the likely fiux based on the colors of ULIRGs 
at similar redshifts and the measured fiuxes at 60 /im. 
We chose cycle times intended to yield signal-to-noise 
ratios of S/N = 50 — 100 for the low- resolution modules 
and S/N > 10 for the high-resolution modules, allowing 
for accurate measurement of faint emission and absorp- 
tion features, as well as accurate spectral decomposition. 

For the 24 objects in the dedicated OHM program, we 
took dedicated sample-up-the-ramp (SUR) peakup ob- 
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TABLE 1 

Radio, optical, and FIR properties of OHMs and non-masing galaxies 
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0.1076 


500 


11.39 
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1037 
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^ Heliocentric optical redshift dDarling fc Giovane ll]| |2002an . 

Computed according to the prescription of lSanders & Mirabell (]1996ll . with a scale factor of C — 1.6. IRAS photometry is from ISanders et al.l i2003D : 
a range in Lfir means that the object was not detected by IRAS at ICQ //m. 
^ OH fluxes are from [Pearling & Giovarie ig I l2002ai rbl): limits are computed according to Eqn. fl] 

Flux densities at 1.4 GHz are from the NRAO VLA Sky Survey jGondon et aI.lfT99^ . 
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servations in both the blue (16 fim) and red (22 /im) 
filters in sample-up-the-ramp (SUR) mode for photomet- 
ric calibration of the spectra. The SUR peakups did not 
coincide in time with the spectroscopy for the major- 
ity of targets. None of the galaxies selected from the 
Spitzer archive possessed SUR peakup data; 16 archived 
galaxies have double correlated sampling (DCS) peakups 
with slightly worse photometric accuracy than those with 
SUR data. Only one galaxy, IRAS 10173-H0828, had 
no peakup observations of any kind. Measured peakup 
fluxes are given in Table [31 

The Spitzer beam is diffraction limited past 6 /xm, with 
slit widths for the spectral modules between 3 — 11". The 
vast majority of the nuclei in merging ULIRGs have an- 
gular separations less than the instrument PSF, and thus 
can be treated as point sources. For the few galaxies close 
enough to be resolved, we chose staring mode observa- 
tions centered on the IR-dominant nuclear region of the 
galaxy. The only available observations of the double 
nucleus in IRAS 10485—1447 were centered on the west- 
ern nucleus. Details of the IRS observations are given in 
Table [2 

4. DATA REDUCTION 
4.1. Dedicated observations of OHM galaxies 

The data were processed using the Spitzer Science Cen- 
ter S17.0 data pipeline. We used basic calibrated data 
(BCD) products for our analysis, having already been 
corrected for flat-fielding, stray light contributions, non- 
linear responsivity in the pixels, and "drooping" (an in- 
crease in detector pixel voltage that occurs during non- 
destructive readouts). The 2-D BCD images were first 
medianed over the data cycles at each nod position to 
remove transient effects such as cosmic rays. For the SL 
and LL modules, we subtracted the sky contribution by 
differencing the BCD images for each nod position with 
the adjacent position in the same module. 

The slit sizes of the SH and LH modules are too 
small to permit extraction of a sky background dur- 
ing the same observation. The continuum levels in the 
high-resolution modules, however, contain strong contri- 
butions from scattered zodiacal light. Estimations of 
the flux at 15 /.tm using SPOT predict zodi contribu- 
tions within the Spitzer beam ranging from 20-80 mJy, 
which in many cases is of comparable magnitude to 
the expected signal from the galaxies themselves. The 
wavelength-dependent brightness of the sky contribution 
means that it cannot be corrected using a simple scaling, 
and so we did not attempt to further calibrate either HR 
module. The calibrated low-resolution spectra were thus 
used for absolute fluxes and the high-resolution spectra 
for line ratio diagnostics (which are unaffected by con- 
tinuum levels). 

To obtain more accurate measurements of faint lines at 
longer wavelengths, we took dedicated off-source sky ob- 
servations in the LH module for all 24 OHMs in our pro- 
gram. Subtraction of the wavelength-dependent back- 
ground, however, significantly affected the measured line 
fluxes. For galaxies with background subtraction in only 
the LH module, this changes the line ratios measured be- 
tween different HR modules [eg, [S III]A33/[S IV]A10.5). 
We reduced the data both with and without subtraction 
of the sky backgrounds; the background-subtracted LH 



spectra are attached as an appendix. 

Following the initial cleaning of the 2-D BCD prod- 
ucts, we eliminated rogue pixels using the IDL package 
IRSCLEAN_MASKEI. We used rogue pixel masks pro- 
vided by the SSC for each IRS campaign, and supple- 
mented the standard masks with manual cleaning of each 
nod and module. The 1-D spectra were then extracted 
using the Spitzer IRS Custom Extractor (SPICE) v. 2.0. 
For all modules we used the optimal extraction routine 
with the standard aperture to improve the S/N ratio in 
faint galaxies. 

The low-resolution modules were stitched together to 
match continuum levels by using a multiplicative scal- 
ing. We fixed the LLl module and then scaled LL2 to 
LLl, SLl to LL2, and SL2 to SLl. The mean scaling 
factors were 1.03 ± 0.08 for LL2 to LLl, 1.49 ± 0.69 for 
SLl to LL2, and 1.37±0.78 for SL2 to SLl. We then cal- 
ibrated the entire low-resolution spectra as a single unit 
by scahng to the IRS 22 ^m sample-up-the-ramp (SUR) 
peakups. The required scaling in the majority of cases 
was quite small, indicating that the sky subtraction and 
spectral extraction techniques are robust; the mean scal- 
ing factor was 0.94 ± 0.08. The acc uracy of the overa ll 
continuum fiux calibration is ^ 5% ( Houck et al.l 1200 4!). 

Noisy areas on both ends of the SH and LH orders 
were trimmed from the 1-D spectra. These areas typi- 
cally encompass a range of 10-30 pixels on the edges of 
the orders and correspond to areas of decreased sensitiv- 
ity on the detector. We deliberately trimmed only pixels 
with an overlapping wavelength range in adjacent orders 
so that a maximum amount of information is preserved. 
In isolated cases, we also removed obvious rogue pixels 
by hand from spectra in which exceptional 1-channel fea- 
tures appear in only a single nod. 

We calculated a simple figure of merit to measure the 
signal-to- noise ratio in the low-resolution data. The data 
near Xrest = 21 yum (a feature- free area near the center 
of most spectra) are fit with a low-order polynomial; the 
median fiux in that region is then divided by the rms 
noise to yield the S/N (Table [3|). We note that this pa- 
rameter is a function of wavelength, as well as the per- 
formance and integration time in each spectral module; 
this is intended to give only a rough estimate for each 
object. The S/N for the samples ranges from '^10— 110, 
with a median of 35. 

4.2. Archival data 

Since the archival OHM and non-masing galaxies did 
not come from a unified observing program, the version 
of the Spitzer data pipeline and the level of process- 
ing varied slightly from object to object - we used the 
most recent versions available in the archive (vl5.3.0 or 
later). The reduction process was identical to that for 
the OHM galaxies in our program, with the exception 
of the LR photometric scaling; since observations in the 
archive varied in availability of peakup data, we used a 
variety of sources to calibrate the spectra. In order of pri- 
ority, we used the IRS dedicated 22 /xm SUR peakups, 
IRS acquisition 16 and 22 ^m DCS peakups, or IRAS 
25 fim observations (all lying within the coverage of the 

1 IRSCLEAN.MASK is available at 

http : //ssc . spitzer . caltech . edu/archanaly/ 
contributed/irsclean/IRSCLEAN.MASK.html 
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TABLE 2 

IRS OBSERVATION LOG 



Object Date Peakup type SLl SL2 LLl LL2 SH LH Program 



IRAS 


01355—1814 




BPU-offsct 


60 


X 


2 


60 


X 


2 


30 


X 


2 


30 


X 


2 


120 


X 


4 


240 


X 


5 


4,6 


IRAS 


01418+1651 




BPU 


14 


X 


2 


14 


X 


2 


14 


X 


4 


14 


X 


4 


120 


X 


5 


60 


X 


7 


2,3 


IRAS 


015Di+i5i8 


2007 bep 09 


BPU-offsct 


14 


X 


7 


14 


X 


7 


30 


X 


3 


30 


X 


3 


120 


X 


2 


60 


X 


2 


1 


IRAS 


02524+2046 


<^c\f\^ O r\c\ 

2007 bep 09 


BPU 


14 


X 


7 


14 


X 


7 


30 


X 


3 


30 


X 


3 


120 


X 


2 


60 


X 


2 


1 


IRAS 


03521+0028 


BPU- offset 


60 


X 


2 


60 


X 


2 


30 


X 


3 


30 


X 


3 


120 


X 


3 


60 


X 


4 


4 


iHAo 


U41zi-t-U2zo 


20U7 Uct Uo 


hSFU-oiiset 


60 


X 


2 


60 


X 


2 


30 


X 


3 


30 


X 


3 


120 


X 


2 


60 


X 


2 


1 


TO A C 

iHAb 


04454— 4ooo 




BPU 


60 


X 


2 


60 


X 


2 


30 


X 


3 


30 


X 


3 


120 


X 


3 


240 


X 


2 


3 


IRAS 


06487+2208 


2007 May 04 


BPU -onset 


14 


X 


7 


14 


X 


7 


30 


X 


3 


30 


X 


3 


120 


X 


2 


60 


X 


2 


1 


IRAS 


07163+0817 


2007 May 03 


BPU 


14 


X 


6 


14 


X 


6 


30 


X 


3 


30 


X 


3 


120 


X 


2 


60 


X 


2 


1 


IRAS 


07572+0533 


2007 May 04 


BPU 


14 


X 


6 


14 


X 


6 


30 


X 


3 


30 


X 


3 


120 


X 


2 


60 


X 


2 


1 


IRAS 


08201+2801 


or\r»T TV /f , no 

2007 May 03 


BPU -onset 


11 


X 


7 


14 


X 


7 


30 


X 


3 


30 


X 


3 


120 


X 


2 


60 


X 


2 


1 


IRAS 


08449+2332 


2007 May 04 


TDT)TT . . IT. . - . 1 

BPU -onset 


14 


X 


6 


14 


X 


6 


30 


X 


3 


30 


X 


3 


120 


X 


2 


60 


X 


2 


1 


IRAS 


r\0 A'7 A 1 1 O 1 o 

08474+1813 


2007 Dec 05 


BPU 


14 


X 


6 


14 


X 


6 


30 


X 


3 


30 


X 


3 


120 


X 


2 


60 


X 


2 


1 


IRAS 


09039+0503 




BPU 


60 


X 


2 


60 


X 


2 


30 


X 


4 


30 


X 


4 


120 


X 


4 


240 


X 


3 


5,6 


IRAS 


09539+0857 




BPU 


60 


X 


2 


60 


X 


7 


30 


X 


4 


30 


X 


4 


120 


X 


3 


240 


X 


3 


5,6 


IRAS 


10035+2740 


onnr^ T nn 

2007 Jun 09 


BPU -onset 


14 


X 


7 


14 


X 


7 


30 


X 


3 


30 


X 


3 


120 


X 


2 


60 


X 


2 


1 


IRAS 


10039—3338 




BPU 


14 


X 


6 


14 


X 


6 


14 


X 


4 


14 


X 


4 


30 


X 


6 


60 


X 


2 


3 


IRAS 


10173+0828 




none 


60 


X 


1 


60 


X 


1 


14 


X 


4 


14 


X 


4 


120 


X 


4 


60 X 12 


3,7,8 


IRAS 


10339+1548 


2007 Jun 08 


BF U-oiiset 


14 


X 


7 


14 


X 


7 


30 


X 


3 


30 


X 


3 


120 


X 


2 


fin 
DU 


X 


o 
z 


1 


IRAS 


10378+1109 




BPU-oiiset 


60 


X 


2 


60 


X 


2 


30 


X 


3 


30 


X 


3 


120 


X 


3 


60 


X 


4 


4 


IRAS 


A o A A^ 

10485—1447 




BPU 


60 


X 


2 


60 


X 


2 


30 


X 


4 


30 


X 


4 


120 


X 


2 


240 


X 


2 


5,6 


TO A C 


ll[j2o-\-61ov 


2007 Jun 09 


BPU -onset 


60 


X 


2 


60 


X 


2 


30 


X 


3 


30 


X 


3 


120 


X 


2 


60 


X 


2 


1 


TO A C 

IRAd 


111 on 1 1 f^oQ 


2007 Jun Oo 


BJr^U -onset 


1 /I 

14 


X 


n 

I 


14 


X 


■7 


30 


X 


3 


30 


X 


3 


120 


X 


2 


60 


X 


2 


1 


IRAS 


ii524+i05o 


onrvT T 1 o 

2007 Jun 12 


BPU -onset 


14 


X 


7 


14 


X 


7 


30 


X 


3 


30 


X 


3 


120 


X 


2 


60 


X 


2 


1 


IRAS 


12018+1941 




BPU- offset 


60 


X 


1 


60 


X 


1 


30 


X 


3 


30 


X 


3 


120 


X 


3 


60 


X 


4 


4 


IRAS 


12032+1707 




BPU- offset 


60 


X 


2 


60 


X 


2 


30 


X 


2 


30 


X 


2 


120 


X 


3 


240 


X 


3 


4,6 


IRAS 


12112+0305 




BPU- offset 


14 


X 


3 


14 


X 


3 


30 


X 


2 


30 


X 


2 


120 


X 


2 


60 


X 


4 


4 


IRAS 


12540+5708 




BPU 


14 


X 


2 


14 


X 


2 


6 


X 


5 


6 


X 


5 


30 


X 


6 


60 


X 


4 


9 


IRAS 


13218+0552 




BPU -onset 


60 


X 


1 


60 


X 


1 


30 


X 


3 


30 


X 


3 


120 


X 


3 


60 


X 


4 


4 


IRAS 


13428+5608 




BPU 


14 


X 


2 


14 


X 


2 


14 


X 


2 


14 


X 


2 


30 


X 


6 


60 


X 


4 


4 


IRAS 


13451+1232 




BPU 


14 


X 


3 


14 


X 


3 


30 


X 


2 


30 


X 


2 


30 


X 


6 


60 


X 


4 


4 


IRAS 


14059+2000 


2007 Jul 31 


BF U-oitset 


60 


X 


2 


60 


X 


2 


30 


X 


3 


30 


X 


3 


120 


X 


2 


60 


X 


2 


1 


TO A C 


14070+0525 




BF U-oliset 


60 


X 


2 


60 


X 


2 


30 


X 


2 


30 


X 


2 


120 


X 


3 


240 


X 


2 


4 


IRAS 


14553+1245 


2007 Jul 31 


BPU 


60 


X 


2 


60 


X 


2 


30 


X 


3 


30 


X 


3 


120 


X 


2 


60 


X 


2 


1 


IRAS 


15327+2340 




BPU 


14 


X 


3 


14 


X 


3 


6 


X 


5 


6 


X 


5 


30 


X 


6 


60 


X 


4 


10 


IRAS 


16090—0139 




T~> TIT T CC i 

BFU-ofiset 


60 


X 


1 


60 


X 


1 


30 


X 


3 


30 


X 


3 


120 


X 


2 


60 


X 


4 


4 


IRAS 




onOR e'en-* 1 7 


BPU 


fin 


^ 


2 


fin 


^ 


2 


OU 


^ 


Q 

o 


OU 


^ 


Q 

o 


120 


^ 


2 


60 


X 


2 


]^ 




iDoUU+ioof5 




rjr^ U-oiiset 


fin 


X 


o 
z 


fin 
6U 


X 


z 


Qn 
3U 


X 


tr 



Qn 
3U 


X 


c: 



1 on 
IzU 


X 


A 

4 


240 


X 


4 


/I fi 

4,0 


TO A Q 

IRAb 


1 ^onv c\c\ 1 /I 
1 /2U i — 0U14 




OOTT 


14 


X 


3 


14 


X 


3 


14 


X 


3 


30 


X 


2 


30 


X 


6 


60 


X 


4 


4 


TO A Q 


1 QQAQ 1 Q /I n 

iooDo-|-oo4y 


zuu ( May Ui 


OOT T 


fin 


X 


z 


fin 
6U 


X 


z 


Qn 
3U 


X 


3 


Qn 
3U 


X 


Q 

3 


1 on 
IzU 


X 


o 
z 


60 


X 


2 


1 


TO A C 


1 Q C^QQ 1 Qf^ 1 7 

ioOoo+ooi / 


zUUD iNOV ZU 


OOT T 


fin 


X 


o 
z 


fin 
60 


X 


z 


Qn 
3U 


X 


Q 

3 


Qn 
3U 


X 


Q 

3 


1 on 
IzU 


X 


o 
z 


60 


X 


2 


1 


TO A Q 

IRAb 


on 1 nn /i i c^j^ 
2U1U0— 4iOD 




OOTT .^fT'r^^+ 

rSFU-ottset 


60 


X 


1 


60 


X 


1 


30 


X 


2 


30 


X 


2 


120 


X 


2 


60 


X 


4 


4 


TO A Q 

IRAb 


2U2oD+1o4d 


onn^j M-->, T on 
200d iNOV ZO 


OOT T 


60 


X 


2 


60 


X 


2 


30 


X 


3 


30 


X 


3 


120 


X 


2 


60 


X 


2 


1 


TO A C 

IKAb 


210/ ( +0000 


onnT T. .n 1 Q 
200; Jun io 


OOT T 

rSFU 


60 


X 


2 


60 


X 


2 


30 


X 


3 


30 


X 


3 


120 


X 


2 


60 


X 


2 


1 


TO A C 

iHAb 


OI OTO 1 Od y1 
21Z / Z+Z0i4 




OOTT ^flPn«4- 

BFU-oltset 


60 


X 


2 


60 


X 


2 


30 


X 


3 


30 


X 


3 


120 


X 


2 


240 


X 


1 


4,11 


TO A C 

iHAb 


2z0oo+oU24 


OnA'7 Tni-i OT 

2007 Jun 27 


OOT T 

BFU 


60 


X 


2 


60 


X 


2 


30 


X 


3 


30 


X 


3 


120 


X 


2 


60 


X 


2 


1 


TO A Q 
ixlAb 


ZZ±iD+U4o ( 


ZUUD JJec Zi 


T30TT ^flPci^4- 

BF u-ortset 


fin 
DU 


X 


z 


fin 
6U 


X 


o 
z 


Qn 
3U 


X 


Q 
O 


Qn 
3U 


X 


Q 

3 


1 on 
IzU 


X 


o 
z 


60 


X 


2 


i 


TO A C 

ixlAb 


Zz4y ± — ±oUo 




BF u-ortset 


fin 
DU 


X 


1 
i 


fin 
6U 


X 


1 


Qn 
3U 


X 


z 


Qn 
oU 


X 


O 

z 


1 on 
IzU 


X 


o 
z 


60 


X 


4 


4 


TO A C 
IRAb 


2aU2o+U72o 


200d Uec 20 


OOTT 

BFU 


1 A 

14 


X 


7 


1 A 

14 


X 


1 


Qn 
3U 


X 


3 


Qn 
3U 


X 


Q 

3 


120 


X 


2 


60 


X 


2 


1 


TO A G 
IKAb 


Zozoo+Uy4D 




OOT T 


fin 
DU 


X 


z 


fin 
6U 


X 


o 
z 


Qn 
3U 


X 


4 


Qn 
3U 


X 


A 

4 


1 on 
IzU 


X 


t) 


240 


X 


4 


fi 
0,6 


TO A G 
IrxAb 


ZOODO + ODU4 




OOT T 


14 


X 


3 


1 A 

14 


X 


3 


Qn 
3U 


X 


z 


Qn 
3U 


X 


z 


Qn 
3U 


X 


fi 
6 


60 


X 


4 


A 

4 


IRAb 


00163—1039 




— TTrm 

BFU 


14 


X 


3 


14 


X 


3 


14 


X 


2 


14 


X 


2 


30 


X 


3 


60 


X 


2 


3 


IRAS 


01572+0009 




BPU 


14 


X 


3 


14 


X 


3 


30 


X 


2 


30 


X 


2 


30 


X 


6 


60 


X 


4 


4 


IRAS 


05083+7936 




BPU 


14 


X 


6 


14 


X 


6 


14 


X 


4 


14 


X 


4 


30 


X 


6 


60 


X 


2 


3 


IRAS 


06538+4628 




BPU 


14 


X 


3 


14 


X 


3 


14 


X 


2 


14 


X 


2 


30 


X 


3 


60 


X 


2 


3 


IRAS 


08559+1053 




BPU-offset 


60 


X 


2 


60 


X 


2 


30 


X 


3 


30 


X 


3 


120 


X 


3 


60 


X 


4 


12 


IRAS 


09437+0317 




BPU 


60 


X 


2 


60 


X 


2 


30 


X 


3 


30 


X 


3 


120 


X 


3 


240 


X 


2 


3 


IRAS 


10565+2448 




BPU 


14 


X 


3 


14 


X 


3 


30 


X 


2 


30 


X 


2 


30 


X 


6 


60 


X 


4 


4 


IRAS 


11119+3257 




BPU 


60 


X 


1 


60 


X 


1 


30 


X 


3 


30 


X 


3 


120 


X 


3 


60 


X 


4 


4 


IRAS 


13349+2438 




BPU 


14 


X 


5 


14 


X 


5 


14 


X 


5 


14 


X 


5 


120 


X 


5 


60 X 10 


13 


IRAS 


15001+1433 




BPU-oflfset 


60 


X 


2 


60 


X 


2 


30 


X 


3 


30 


X 


3 


120 


X 


3 


60 


X 


4 


4 


IRAS 


15206+3342 




BPU-offset 


60 


X 


1 


60 


X 


1 


30 


X 


3 


30 


X 


3 


120 


X 


3 


60 


X 


4 


4 


IRAS 


20460+1925 




BPU 


14 


X 


5 


14 


X 


5 


14 


X 


5 


14 


X 


5 


120 


X 


5 


60 X 10 


13 


IRAS 


23007+0836 




BPU 


14 


X 


2 


14 


X 


2 


6 


X 


5 


6 


X 


5 


30 


X 


4 


60 


X 


2 


14 


IRAS 


23394-0353 




BPU 


30 


X 


6 


30 


X 


6 


60 


X 


2 


60 


X 


2 


14 


X 


6 


14 


X 


4 


3 


IRAS 


23498+2423 




BPU-offset 


60 


X 


2 


60 


X 


2 


30 


X 


2 


30 


X 


2 


120 


X 


3 


240 


X 


2 


1 



Note. — Spitzer archival data are from programs: (1) 30407 (PI: J. Darling); (2) - 3237 (PI: E. Sturm); (3) - 30323 (PI: 
L. Armus); (4) - 105 (PI: J. Houck); (5) - 2306 (PI: M. Imanishi); (6) - 3187 (PI: S. Veilleux); (7) - 3605 (PI: C. Bradford); (8) 
- 20549 (PI: R. Joseph); (9) - 1442 (PI: L. Armus); (10) - 1444 (PI: L. Armus); (11) - 20375 (PI: L. Armus); (12) - 666 (PI: 
J. Houck): (13) - 61 (PI: G. Rieke); (14) - 14 (PI: J. Houck). Exposure times for all modules are given as seconds per cycle X 
number of cycles. 
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LL modules). For galaxies from the archive, 17 are cali- 
brated with DCS peakups, 21 with IRAS 25 /zm photom- 
etry, and 4 are left uncalibrated. The mean scaling factor 
for the objects without SUR photometry was 1.05 ±0.25, 
slightly higher than the mean scaling for galaxies in the 
dedicated OHM program. 

Five OHMs and five non-masing galaxies from the 
archive had both SH and LH sky backgrounds taken 
simultaneously with the spectroscopic observations; the 
remainder had no high-resolution sky backgrounds in ei- 
ther module. Since we emphasize uniformity of the ob- 
servations to the fullest extent possible, all data used 
for statistical comparisons between the samples use data 
without HR sky subtraction; line measurements for the 
background-subtracted objects are given in the Ap- 
pendix. 

The spectra for IRAS 20460+1925 and 
IRAS 23028-^0725 had no flux in the SL modules, 
most likely due to a pointing error during observations. 
No SL data from either galaxy are used in our analyses 
here or in Paper II. 

5. RESULTS 

We show examples of the peakup-scaled low-resolution 
spectra for the OHMs in Figure [H with the individual 
modules stitched together and bonus orders removed. 
Examples of the high-resolution SH and LH data are 
shown in Figures [3] and HJ full spectra for all galaxies 
are available as an online supplement. While individ- 
ual orders within the high-resolution modules are typ- 
ically well-aligned in flux, the differences in calibration 
between the SH and LH modules are clearly apparent 
when matching the spectra; this is due to a combination 
of different slit sizes for the SH and LH modules (a fac- 
tor of ~ 4) and the lack of separate sky subtraction for 
the SH modules. For this reason, as well as emphasizing 
the narrow atomic and molecular features visible in the 
high-resolution spectra, we display separate plots for the 
SH and LH modules. 

A small amount of the archival objects have pre- 
viously p ublished full IRS spectra (lArmus et al.l 120041 : 
IWeedman et all l2005t lArmus et all I2007D. maiii l y con - 
sisting of bright, nearby galaxies. iFarrah et al.l ()2007f ) 
published HR spectra for roughly half of our archival 
OHM hosts. While many papers use the available data 
from the GTO programs, however, the majority of ob- 
jects extracted from the archive have no published spec- 
tra, although some da ta are used in larg e r studies of 
ULIR G properties (eg , iHigdon et al.l 120061 : iDesai et"al] 
120071: iHao et al.ll2007r) . The spectra for many of the 
archival galaxies are thus presented here for the first 
time. 

Comparison of our data with the few spectra of objects 
previously published {eg, Mrk 1014, Arp 220, NGC 7469) 
revealed no significant differences in spectral shape or 
detection of individual features. Measurements of line 
flux and equivalent widths, however, may be affected by 
the photometric scaling and/or line-fitting routines used; 
for this reason, we chose to reduce all data in a uniform 
matter. 

For many of our high-resolution spectra, especially 
those with low S/N, there exist individual spikes 
that do not correspond to any identified feature (see 
IRAS 01562-f 2528 for a prominent example). These fea- 



tures are typically 1-2 channels wide, much narrower 
than the expected linewidth for an unresolved feature. 
We regard such features as spurious, possibly caused by 
hot pixels or other instrumental conditions that are not 
corrected by our cleaning routines. All features we regard 
as valid detections are listed in the data tables, with the 
locations of the most common features marked on the 
spectra themselves. 

5.1. Continuum 

The continuum emission for all objects in the OHM 
sample has a relatively homogenous spectral shape over 
the range of the IRS, although differences in spectral 
shape between the two samples do appear and are ex- 
plored in Paper II. Figure [5] shows the individual ob- 
jects overlaid with a template generated by median- 
ing the flux in each wavelength bin from all galaxies. 
The template bears a close resemblanc e to starbursting 
ULIRG spectra seen in pre vious surveys (|Hao et al.ll2007l : 
IWeedman fc Houc3 120091) . The LR spectrum clearly 
shows silicate absorption at 9.7 and 18 /xm and water 
ice absorption at 6 fim. Low-resolution emission features 
are dominated by the broad PAH features from 6-13 /im, 
with weaker contibutions from neon, sulfur, and molec- 
ular hydrogen also visible. 

The continuum data from Xrest — 20 — 30 fim is in 
most cases well-characterized by a power-law fit, with 
the short wavelength break occurring near the 18 /xm sil- 
icate feature and the long wavelength end cut off by the 
spectral range of the IRS. Shortward of 15 /im, the con- 
tinuum becomes increasingly contaminated by individual 
absorption and emission features, especially from PAH 
emissio n and the dee p silica te absorption at 9.7 /im. Fol- 
lowing iBrandl et all (|2006[ ). we fit a spectral index to 
the continuum in two components, with Q!3o_2o measur- 
ing the relatively feature-free flux from 20 — 30 /im and 
q;i5_6 measuring the contribution from 5.3 — 14.8 /im; the 
wavelengths are slightly shifted to avoid contamination 
from water ice at 6 /im and [Ne III] at 15.6 /im (Table|3|). 

The mean 15-6 spectral index for the entire OHM 
sample is q;i5-6 = 2.1 ± 0.6; the mean 30-20 slope 
is q;3o_2o — 4.8 ± 1.1. The shallowest 15-6 slope oc- 
curs for IRAS 10039-3338 (ais^e = -0.1), an object 
with weak PAH features and very strong silicate absorp- 
tion; the steepest index occurs for IRAS 12018+1941 
(q;i5_6 — 3.0), which has moderate PAH and line emis- 
sion features and a nearly constant spectral index over 
the entire mid-IR range. Steeper 15-6 indices are likely 
due to a combination of smaller relative quantities of 
warm dust (thermal blackbodies of ^ 300 K peaking near 
10 /im) and larger quantities of cooler dust. 

The shallowest 30-20 slope occurs for 
IRAS 13451+1232 (030-20 = 1-9); the low-resolution 
spectrum for this object more closely re sembles that seen 
in Seyfert galaxies an d PG quasars ([Schweitzer et al.l 
120061: IHao et all I2007D . with weak PAH emission and 
shallower silicate absorption. The continuum emission 
for this object is also much closer to being uniform over 
the entire range of the IRS; the difference between the 
two spectral indices is only Aa = —0.4, compared to an 
average of Aa = 2.7 for the entire sample. This behavior 
is more typical of non-thermal emission that can extend 
over many decades with the same index. The steepest 
30-20 emission measured is from IRAS 11028+3130 
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Fig. 2. — IRS spectra from the low-resolution modules (LR) for OHMs. Spectra for all OHMs and non-masing galaxies are available as 
an online supplement; portions are shown here for guidance on form and style. All detected PAH emission and absorption features from 
water ice and silicates are marked. 
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TABLE 3 

IRS PHOTOMETRY AND CONTINUUM MEASUREMENTS 



Object 16 nm PU 22 /im PU PU type «i5-(i «3o-20 S/N 







[m Tvl 


\m Tvl 










IRAS 


01355-1814 




45.5 


DCS 


2.8 + 0.3 


5.2 + 0.8 


75 


IRAS 


01418+1651 








2.2 ± 0.3 


5.2 + 0.8 


18 


IRAS 


01562+2528 


6.1 


13.0 


SUR 


1.8 + 0.3 


5.1 + 0.8 


11 


IRAS 


02524+2046 


10.6 


21.7 


SUR 


2.0 + 0.4 


5.7+1.0 


25 


IRAS 


03521+0028 


25.6 




DCS 


2.3 + 0.4 


5.5 + 0.8 


63 


IRAS 


04121+0223 


10.5 


22.6 


SUR 


2.6 ± 0.1 


5.9 ± 0.5 


21 


IRAS 


04454-4838 








2.7 + 0.4 


6.1 + 0.9 


37 


IRAS 


06487+2208 


85.0 


177.5 


SUR 


2.4 + 0.2 


3.6 + 0.8 


48 


IRAS 


07163+0817 


12.3 


28.8 


SUR 


2.5 + 0.3 


4.8 + 1.0 


26 


IRAS 


07572+0533 


52.1 


103.0 


SUR 


2.6 + 0.3 


3.0 + 0.7 


55 


IRAS 


08201+2801 


61.8 


74.3 


SUR 


2.2 + 0.2 


4.5 + 0.7 


37 


IRAS 


08449+2332 


26.6 


48.2 


SUR 


2.2 ± 0.2 


4.4 + 0.4 


52 


IRAS 


08474+1813 


9.5 


28.7 


SUR 


2.2 ± 0.5 


6.5 + 0.6 


30 


IRAS 


09039+0503 


86.7 




DCS 


1.8 ± 0.5 


5.5 + 0.7 


41 


IRAS 


09539+0857 








2.2 ± 0.1 


6.0 + 0.3 


36 


IRAS 


10035+2740 


9.2 


22.8 


SUR 


2.1 ± 0.2 


6.0 + 0.5 


23 


IRAS 


10039-3338 








-0.1 ± 0.1 


5.4 + 0.3 


17 


IRAS 


10173+0828 








2.6 + 0.4 


6.4 + 1.0 


48 


IRAS 


10339+1548 


10.7 


28.1 


SUR 


2.5 + 0.3 


5.1 + 0.7 


22 


IRAS 


10378+1109 




114.8 


DCS 


2.0 + 0.3 


4.5 + 0.7 


35 


IRAS 


10485-1447 








2.5 + 0.3 


4.8 + 0.8 


75 


IRAS 


11028+3130 


4.8 


13.0 


SUR 


2.4 + 0.2 


7.0 + 0.6 


26 


IRAS 


11180+1623 


14.8 


33.6 


SUR 


2.2 + 0.3 


5.6 + 0.9 


15 


IRAS 


11524+1058 


8.2 


14.8 


SUR 


1.5 + 0.4 


6.1 + 0.8 


20 


IRAS 


12018+1941 


121.7 




DCS 


3.0 + 0.3 


2.9 + 0.4 


63 


IRAS 


12032+1707 




73.9 


DCS 


2.2 ± 0.3 


4.2 + 0.7 


42 


IRAS 


12112+0305 


91.9 




DCS 


2.3 ± 0.3 


5.3 + 0.7 


42 


IRAS 


12540+5708 








1.6 ± 0.3 


2.4 + 0.7 


30 


IRAS 


13218+0552 


212.9 




DCS 


0.6 ± 0.3 


2.4 + 0.7 


100 


IRAS 


13428+5608 








1.9 ± 1.3 


4.6 + 0.6 


12 


IRAS 


13451+1232 








2.3 ± 0.4 


1.9 + 0.7 


50 


IRAS 


14059+2000 


12.2 


26.3 


SUR 


1.1 ± 0.6 


4.9 + 0.8 


20 


IRAS 


14070+0525 




30.1 


DCS 


1.6 + 0.3 


5.6 + 0.5 


33 


IRAS 


14553+1245 


28.7 


61.2 


SUR 


2.0 + 0.7 


4.3 + 0.7 


22 


IRAS 


15327+2340 








2.4 + 0.2 


5.5 + 0.4 


17 


IRAS 


16090-0139 


71.4 




DCS 


1.3 + 0.4 


4.6 + 0.6 


59 


IRAS 


16255+2801 


16.0 


36.7 


SUR 


1.4 + 0.3 


4.7 + 0.6 


22 


IRAS 


16300+1558 




35.0 


DCS 


1.7 + 0.5 


5.6 + 0.9 


31 


IRAS 


17207-0014 








1.8 ± 0.6 


6.2 + 0.6 


28 


IRAS 


18368+3549 


26.6 


49.8 


SUR 


1.7 + 0.5 


5.4 + 0.8 


29 


IRAS 


18588+3517 


43.6 


92.1 


SUR 


1.5 ± 0.3 


4.8 + 0.8 


28 


IRAS 


20100-4156 


86.7 




DCS 


1.9 ± 0.2 


5.2 + 0.7 


73 


IRAS 


20286+1846 


11.4 


22.0 


SUR 


2.3 ± 0.3 


5.8 + 0.6 


20 


IRAS 


21077+3358 


31.1 


62.1 


SUR 


2.8 ± 0.3 


4.2 + 0.6 


52 


IRAS 


21272+2514 




39.0 


DCS 


2.0 + 0.3 


5.0 + 0.7 


36 


IRAS 


22055+3024 


54.5 


132.4 


SUR 


2.7 + 0.3 


4.0 + 0.7 


30 


IRAS 


22116+0437 


46.6 


68.2 


SUR 


2.3 + 0.5 


4.5 + 0.9 


35 


IRAS 


22491-1808 


87.5 




DCS 


2.8 + 0.3 


5.1 + 0.6 


28 


IRAS 


23028+0725 


56.4 


140.7 


SUR 




3.2 + 0.7 


25 


IRAS 


23233+0946 








2.1 + 0.3 


4.7 + 0.5 


33 


IRAS 


23365+3604 








2.5 + 3.1 


4.6 + 0.5 


15 


IRAS 


00163-1039 








2.3 + 0.4 


2.1 ±0.4 


34 


IRAS 


01572+0009 








1.8 + 0.2 


2.2 + 0.3 


81 


IRAS 


05083+7936 








2.0 + 0.4 


2.8 + 0.4 


48 


IRAS 


06538+4628 








2.9 ± 0.3 


2.5 + 0.4 


37 


IRAS 


08559+1053 




90.7 


DCS 


1.3 + 0.1 


2.8 + 0.4 


64 


IRAS 


09437+0317 








1.7 + 0.3 


2.2 + 0.4 


19 


IRAS 


10565+2448 








2.1 ± 0.4 


3.2 + 0.5 


19 


IRAS 


11119+3257 








1.2 ± 0.1 


2.5 ± 0.3 


79 


IRAS 


13349+2438 








0.8 ± 0.1 


0.1 + 0.1 


110 


IRAS 


15001+1433 




135.3 


DCS 


1.9 + 0.2 


3.4 + 0.5 


58 


IRAS 


15206+3342 


110.7 




DCS 


2.3 + 0.4 


2.6 + 0.4 


65 


IRAS 


20460+1925 










1.1 + 0.3 


16 


IRAS 


23007+0836 








1.7 + 0.1 


1.8 + 0.4 


19 


IRAS 


23394-0353 








1.8 + 0.3 


2.9 + 0.4 


37 


IRAS 


23498+2423 


46.5 




DCS 


1.0 + 0.1 


3.1 + 0.4 


87 



Note. — 



Errors in the peakup (PU) fluxes are at the 15% level. 
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Fig. 4. — IRS spectra from the long-high module (LH) for OHMs. Spectra for all OHMs and non-masing galaxies are available as an 
online supplement; portions are shown here for guidance on form and style. All detected atomic and H2 emission features are marked. 
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Fig. 5. — Low-resolution spectrum of all galaxies in our sample 
normalized at S\ = 15 fim. The black spectrum is the composite 
template made from an error-weighted median of the individual 
galaxies; the grey shaded area shows the l-tr envelope for each res- 
olution element (exaggerated toward negative values in log space). 

(a30-20 = 7.0); the galaxy shows moderate PAH and 
hne emission features, but a flat continuum between 12 
and 20 /im. 

The non-masing galaxies have average spectral indices 
of ai5_6 = 1.8 ± 0.6 and a3o_2o = 2.5 ± 0.9. A full sta- 
tistical analysis of the values for the two samples (par- 
ticularly a3o_2o) is presented in Paper II. 

5.2. Atomic emission lines 

We measured emission from atomic and molecular 
lines using the standard packages in the Spectroscopic 
Modehng Analysis a nd Reduction Tool (SMART) v6.2.4 
(jHigdon et al.ll2004l) . A simple Gaussian is a good fit for 
virtually all high-resolution lines in the sample; in cases 
where lines are blended (such as the [Ne V]/[C1 II] and 
[O IV] /[Fe II] complexes), we used a multi-Gaussian fit 
centered at the redshifted rest wavelengths of the ex- 
pected transitions. To compute upper limits for non- 
detections, we use the S-tr noise measured from the 
surrounding continuum and a Gaussian shape with an 
FWHM estimated from detected lines. Accuracy for all 
measured line fiuxes is on the order of ^ 10% (Tables |4] 
andE]). 

The most common lines detected are the forbidden 
[Ne II] A12.814 and [Ne III] A15.555 transitions. [Ne II] is 
observed in nearly the entire sample, with detections 
in 50/51 OHMs and 14/15 non-masing galaxies. The 
only exceptions are the OHM IRAS 11028-1-3130 and the 
non-maser IRAS 20460-1-1925. [Ne III] is also common, 
detected in 43/51 OHMs and 14/15 non-masing galax- 
ies. Other common lines are the [S III] A18.713 (de- 
tected in - 80% of galaxies) and [S IV] A10.511 (- 50%). 
[Ar III] A8.991 is detected in 15 OHMs and 5 non-masing 
galaxies, but the redshifted line is not visible in the SH 
module for archived objects at z < 0.1. 

We detect "rarer" line transitions that appear in less 
than 15% of our sample, including [Ne V] A14.322 and 
A24.318, [Fe II] A17.936 and A25.988, and [O IV] A25.890. 
IRAS I0339-I-1548 is the only galaxy in the sample with 
detections for all atomic transitions listed above. Eight 
galaxies also show the non-forbidden HI 7-6 transition 



(Humphreys-a) at 12.368 fim. 

[S III]A33.481 and [Si II]A34.815 are commonly ob- 
served in ULIRGs for which the redshifted t ransitions ap- 
pear w ithin the red edge of the LH module. iFarrah et al.l 
pOOT*) detected [S III]A33.481 in roughly half of aU 
ULIRGs with z < 0.06, all of which show the stronger 
[S III]A18.713 associated transition. For our sample, 
however, many galaxies have the lines redshifted beyond 
the range of the IRS; when visible, the lines lie at the 
far red edge of the II"^ order in the LH module, an area 
with high noise and decreased sensitivity with respect to 
neighboring orders. 

For galaxies taken from the Spitzer archive, we com- 
pared our high-resoluti on measurements wi th those ap- 
pearing in the sample of lFarrah et al.l (|2007[) . The agree- 
ment between detected lines is go od for nearly all ob- 
jects; however, iFarrah et al.l ()2007l ) report detections of 
[Ne V] in IRAS 11119-^3257, [O IV] in IRAS 13451-1-1232, 
and Ha S(0) and S(2) lines in IRAS 0I572-H0009 and 
IRAS 23498-1-2423 which we fail to confirm. Limits for 
the measured fluxes are given in Tables |4] and [5] and [6] 

The IRS is designed to make accurate measurements 
of narrow atomic transitions in the SH and LH mod- 
ules; however, several lines are also detected in the 
LR modules, most often the powerful neon, sulfur, and 
H2 transitions. The only emission lines visible in the low- 
resolution modules without a corresponding detection in 
high-resolution are the H2 S(7) line at 5.5 fiia and the 
H2 S(5)/[Ar II] complex near 6.7 /im; this is due to the 
SH low-wavelength cutoff at Xrest = 8.2 — 9.0 ^m, de- 
pending on the redshift of the galaxy. All other atomic 
emission features observed in the low-resolution modules 
have corresponding detections in high-resolution; fur- 
thermore, blending of narrow lines makes accurate mea- 
surements of flux difficult in the LR modules. For iso- 
lated lines with well-defined surrounding continuum, our 
fluxes are consistent for measurements in both low- and 
high-resolution. 

We note the presence of two features which have 
no obvious identifications occurring in the LH spec- 
tra for multiple objects: one is an emission feature 
seen near 29 /im (a prominent example occurs for 
IRAS 17539-1-2935) and the second is an absorption fea- 
ture near 30 /im. The two are often paired and are seen in 
^ 50% of the galaxies observed. The rest wavelengths of 
the transitions, however, vary significantly from object 
to object (with a standard deviation of a\ ~ 0.7 /^m), 
while the observed wavelengths are nearly fixed {ax < 
0.05 /xm). This implies that the features are either arti- 
facts of the extraction process or that both the emission 
and absorption come from unidentified foreground fea- 
tures with little to no Doppler shift. Given that we see 
no evidence for these unidentified lines in any of the LR 
spectra (which should be detectable, given the high S/N 
for many of the features), we consider both to be spuri- 
ous. 

5.3. Molecular hydrogen 

We detected multiple emission lines from the pure ro- 
tational series of molecular hydrogen in both OHMs and 
non-masing galaxies. At redshifts of 2: < 0.1, transitions 
from H2 S(0) at 28.22 fj.m to H2 S(3) at 9.67 /zm are 
visible in the HR modules; in addition, the LR module 
is capable of detecting lines as far out as the S(7) tran- 
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TABLE 4 

Atomic line fluxes for high-resolution spectra: 8 — 16 microns 





Object 


[Ar 


III] 


[S 


IV] 


H 


.1 7-6 


— 

[Nc ii] 


[Nc v] 


[CI II] 


— 

[Nc III] 


Xrest [mH 


8.991 


10.511 


12.368 


12.814 


14.322 


14.369 


15.555 


IRAS 


01355—1814 


0, 


.15 




0.10 


< 


0.20 


2.45 


< 


0.27 


< 


0.25 


0.89 


IRAS 


01418+1651 








0.15 




0.13 


4.17 


< 


0.46 


< 


0.29 


0.41 


IRAS 


01562+2528 


< 0, 


.77 


< 


1.93 


< 


0.30 


1.07 


< 


0.91 


< 


0.75 


0.77 


IRAS 


02524+2046 


< 0. 


.63 


< 


1.01 


< 


0.36 


2.39 


< 


0.68 


< 


0.56 


< 0.97 


IRAS 


03521+0028 


< 0. 


.34 


< 


0.40 


< 


0.41 


2.72 


< 


0.36 


< 


0.24 


1.11 


IRAS 


04121+0223 


0, 


.35 


< 


0.84 


< 


0.41 


1.81 


< 


1.01 


< 


0.99 


< 0.87 


IRAS 


04454—4838 








0.47 


< 


0.74 


2.04 


< 


0.75 


< 


0.61 


0.39 


IRAS 


06487+2208 


1, 


,71 




1.78 


< 


0.81 


11.75 


< 


1.44 


< 


1.19 


9.45 


IRAS 


07163+0817 


0, 


,44 


< 


1.05 


< 


0.22 


3.19 


< 


0.46 


< 


0.39 


0.33 


IRAS 


07572+0533 


< 1. 


,69 


< 


0.87 


< 


0.40 


1.74 


< 


0.79 


< 


0.65 


< 0.89 


IRAS 


08201+2801 


< 0, 


,70 


< 


2.74 


< 


0.64 


2.23 


< 


0.60 


< 


0.41 


0.88 


IRAS 


08449+2332 


< 0, 


,94 


< 


0.71 


< 


0.64 


3.60 


< 


2.03 


< 


1.64 


1.57 


IRAS 


08474+1813 


< 0, 


,76 


< 


0.64 


< 


0.36 


1.26 


< 


1.10 


< 


0.89 


< 1.10 


IRAS 


09039+0503 


0, 


,11 


< 


0.27 


< 


0.85 


3.68 


< 


0.46 


< 


0.38 


1.17 


IRAS 


09539+0857 


< 0, 


,63 


< 


0.43 


< 


0.30 


1.25 


< 


0.86 


< 


0.69 


< 1.88 


IRAS 


10035+2740 


< 0, 


,89 


< 


0.75 


< 


0.48 


1.82 


< 


0.42 


< 


0.31 


0.57 


IRAS 


10039—3338 








0.93 


< 


1.25 


16.69 


< 


1.21 


< 


1.10 


3.93 


IRAS 


10173+0828 






< 


0.29 


< 


0.34 


1.71 


< 


0.33 


< 


0.32 


0.46 


IRAS 


10339+1548 


1, 


,00 




0.55 


< 


0.34 


1.25 




1.24 


< 


1.13 


2.03 


IRAS 


10378+1109 


0, 


,34 


< 


0.25 


< 


0.35 


3.95 


< 


0.61 


< 


0.63 


0.68 


IRAS 


10485-1447 


< 1. 


,14 


< 


0.20 


< 


0.23 


1.99 


< 


0.44 


< 


0.48 


0.36 


IRAS 


11028+3130 


< 1. 


,06 


< 


0.70 


< 


0.28 


< 0.63 


< 


0.44 


< 


0.38 


< 0.74 


IRAS 


11180+1623 


< 0, 


,80 


< 


0.70 


< 


0.36 


1.98 


< 


0.50 


< 


0.39 


0.62 


IRAS 


11524+1058 


< 0, 


,78 


< 


0.56 


< 


0.35 


7.82 


< 


0.43 


< 


0.37 


< 0.77 


IRAS 


12018+1941 


< 0, 


,27 


< 


0.28 


< 


0.59 


2.73 


< 


0.52 


< 


0.38 


0.64 


IRAS 


12032+1707 


< 1. 


,33 




0.13 


< 


0.48 


4.98 




1.39: 


< 


1.51 


1.57 


IRAS 


12112+0305 








0.43 


< 


0.88 


13 06 


< 0.76 


< 


0.55 


3.37 


IRAS 


12540+5708 






< 11.45 


< 


4.40 


17.95 


< 8.05 


< 


6.86 


9 36 


IRAS 


13218+0552 


< 0, 


,66 


< 


0.54 


< 


0.46 


0.84 


< 0.71 


< 


1.23 


< 1.24 


IRAS 


13428+5608 








7.69 


< 


2.42 


41.21 


10.62 


< 


5.37 


29.05 


IRAS 


13451+1232 


0, 


.57 




1.57 


< 


0.70 


4.71 




0.71 


< 


0.88 


4.73 


IRAS 


14059+2000 


< 0, 


.64 




0.35 


< 


0.53 


2.83 


< 


0.88 


< 


1.00 


2.65 


IRAS 


14070+0525 


< 0, 


.19 


< 


0.18 


< 


0.27 


1.33 


< 


0.42 


< 


0.33 


2.01 


IRAS 


14553+1245 


0, 


.47 




0.71 


< 


0.25 


3.24 


< 


1.55 


< 


1.29 


3.10 


IRAS 


15327+2340 






< 


0.79 


< 


3.60 


61.13 


< 


6.29 


< 


5.29 


6.89 


IRAS 


16090-0139 


0, 


.49 




0.10 


< 


0.49 


6.72 


< 


0.65 


< 


0.76 


2.26 


IRAS 


16255+2801 


0, 


.18 




0.20 


< 


0.27 


1.73 


< 


0.82 


< 


0.68 


1.18 


IRAS 


16300+1558 


0, 


,13 


< 


0.15 


< 


0.28 


2.27 


< 


0.33 


< 


0.25 


0.42 


IRAS 


17207-0014 








0.40 


< 


2.12 


38.84 


< 


0.93 


< 


0.69 


8.42 


IRAS 


18368+3549 


< 0, 


,42 


< 


0.81 




0.11 


6.91 


< 


0.34 


< 


0.43 


1.07 


IRAS 


18588+3517 


0, 


,57 




0.37 




0.24 


5.12 


< 


0.40 


< 


0.43 


1.87 


IRAS 


20100-4156 


0, 


,33 




0.23 


< 


0.31 


6.73 


< 


0.65 


< 


0.49 


1.71 


IRAS 


20286+1846 


< 0. 


,56 


< 


0.50 


< 


0.27 


1.67 


< 


0.72 


< 


0.65 


0.37 


IRAS 


21077+3358 


< 0, 


.67 


< 


0.63 


< 


0.39 


3.06 


< 


0.45 


< 


0.36 


1.09 


IRAS 


21272+2514 


< 0, 


.18 


< 


0.16 




0.16 


2.29 


< 


0.31 


< 


0.22 


0.35 


IRAS 


22055+3024 


0, 


.12 


< 


0.67 


< 


0.71 


4.55 


< 


0.76 


< 


0.74 


1.05 


IRAS 


22116+0437 


< 0, 


.47 


< 


0.73 


< 


0.37 


2.25 


< 


0.82 


< 


0.68 


1.33 


IRAS 


22491-1808 








0.41 


< 


0.54 


4.88 


< 


0.55 


< 


0.43 


1.70 


IRAS 


23028+0725 


< 0, 


.55 


< 


0.61 


< 


0.39 


1.81 


< 


0.90 


< 


0.72 


1.04 


IRAS 


23233+0946 


0, 


.37 




0.33 


< 


0.43 


4.93 


< 


0.43 


< 


0.38 


1.06 


IRAS 


23365+3604 






< 


0.36 


< 


0.51 


8.51 


< 


0.63 


< 


0.52 


1.12 


IRAS 


00163-1039 








3.40 




0.43 


80.95 


< 1.96 


< 


0.92 


14.53 


IRAS 







.73 




3.06 




0.72 


6.15 




5.51 


< 


3.23 


10.34 


IRAS 


05083+7936 






< 1.19 


< 


1.57 


49.40 




0.60 




0.48 


7.63 


IRAS 


06538+4628 








0.80 




1.09 


47.39 




0.69 


< 


1.16 


6.07 


IRAS 


08559+1053 


< 0, 


.37 




0.56 


< 


0.43 


8.38 




0.51 




0.35 


1.87 


IRAS 


09437+0317 






< 1.04 


< 


0.61 


8.74 


< 0.64 


< 


0.51 


1.22 


IRAS 


10565+2448 






< 0.76 


< 


1.34 


57.60 


< 1.29 




0.67 


7.65 


IRAS 


11119+3257 


< 0. 


.83 




0.31 


< 


0.78 


2.19 


< 0.71 


< 


0.58 


1.89 


IRAS 


13349+2438 


2, 


.76 




1.66 


< 


0.45 


1.43 




0.81 


< 


0.99 


3.50 


IRAS 


15001+1433 


0. 


.39 




0.28 


< 


0.30 


6.61 




1.08 


< 


0.98 


2.62 


IRAS 


15206+3342 


2, 


.01 




3.82 




0.23 


10.96 


< 0.34 


< 


0.41 


19.87 


IRAS 


20460+1925 


< 0, 


.39 


< 0.66 


< 


0.42 


< 0.44 


< 0.60 


< 


0.51 


< 0.69 


IRAS 


23007+0836 








8.67 


< 


1.84 


179.04 




8.36 


< 


4.23 


33.69 


IRAS 


23394-0353 








0.72 




0.57 


46.75 


< 1.45 


< 


1.08 


7.71 


IRAS 


23498+2423 


0, 


.23 




1.27 


< 


0.23 


3.10 




0.92 


< 


0.91 


7.79 


Note 


. — Line fluxes 


1 are 


given 


in 


10-21 


W 




— indicates 


; that the redshifted line 



wavelength lay outside the range of the IRS. 
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TABLE 5 

Atomic line fluxes for high-resolution spectra: 17 — 35 microns 





Object 


— ■ 

[Fc ii] 




[S III] 


[N 


c v] 




[O IV] 


[Fc ii] 




[S III] 




[Si II] 




17.936 


18.713 


24.318 


25.890 


25.988 


33.481 


34.815 


IRAS 


01355—1814 


< 0.51 


< 1.52 


< 


0.78 


< 1.28 


< 


0.82 






IRAS 


01418+1651 


< 0.67 


0.94 


< 


3.84 


< 3.52 


< 


4.05 


< 24.62 


19.43 


IRAS 


01562+2528 


< 0.77 


0.67 


< 


0.73 


< 2.06 


< 


1.29 






IRAS 


02524+2046 




0.47 


< 


1.04 


< 0.96 


< 


0.59 






IRAS 


03521+0028 


1.70 


1.18 


< 


1.02 


< 0.93 


< 


0.59 






IRAS 


04121+0223 


< 0.92 


0.93 


< 


1.06 


0.97 


< 


1.34 






IRAS 


04454—4838 


< 0.27 


< 1.00 


< 


7.75 


< 9.59 


< 


6.29 


0.90 


< 87.49 


IRAS 


06487+2208 


< 0.87 


4.89 


< 


1.36 


0.48 




1.31 






IRAS 


07163+0817 


< 1.04 


0.81 


< 


1.73 


< 0.78 


< 


0.71 






IRAS 


07572+0533 


< 0.61 


< 1.15 


< 


1.37 


< 1.53 


< 


1.03 






IRAS 


08201+2801 


< 0.84 


0.41 


< 


0.50 


< 0.66 


< 


0.50 






IRAS 


08449+2332 


< 0.60 


1.62 


< 


0.73 


< 0.69 


< 


0.40 






IRAS 


08474+1813 


< 1.09 


0.86 


< 


1.27 


< 3.04 


< 


1.86 






IRAS 


09039+0503 


< 0.85 


0.88 


< 


0.55 


< 1.11 




0.48 






IRAS 


09539+0857 


< 1.82 


< 1.02 


< 


1.41 


< 2.05 


< 


1.36 






IRAS 


10035+2740 


< 0.97 


0.78 


< 


0.95 


< 2.11 


< 


1.37 






IRAS 


10039—3338 


< 0.81 


6.57 


< 14.49 


< 11.09 


< 


7.54 


< 10.49 


101.40 


IRAS 


10173+0828 


< 0.42 


< 1.36 


< 1.33 


< 2.42 


< 


1.61 


2.55 


81.70 


IRAS 


10339+1548 


0.49 


1.17 




0.52 


2.85 




0.89 






IRAS 


10378+1109 


< 0.95 


1.82 


< 0.87 


< 1.99 


< 


1.30 






IRAS 


10485-1447 


< 0.87 


< 1.08 


< 0.91 


< 1.08 


< 


0.66 






IRAS 


11028+3130 


< 0.47 


< 0.71 


< 3.07 


< 4.45 


< 


3.14 






IRAS 


11180+1623 


< 0.71 


0.89 


< 0.85 


< 1.00 


< 


0.69 






IRAS 


11524+1058 


< 0.52 


< 1.61 


< 1.49 


< 0.93 


< 


0.70 






IRAS 


12018+1941 


< 1.01 


< 1.32 


< 1.22 


< 1.71 


< 


1.30 






IRAS 


12032+1707 


< 1.22 


1.18 




0.66 


< 1.50 


< 


1.00 






IRAS 


12112+0305 


0.68 


4.32 


< 1.50 


< 7.56 


< 


4.65 


8.77 




IRAS 


12540+5708 


< 4.40 


< 5.05 


< 22.99 


< 18.84 


< 14.43 


< 49.65 




IRAS 


13218+0552 


< 0.78 


< 2.27 


< 1.85 


< 2.93 


< 1.76 






IRAS 


13428+5608 


< 1.32 


16.67 




9.28 


75.13 


< 29.79 


28.14 


< 134.61 


IRAS 


13451+1232 


< 1.19 


1.15 


< 2.50 


< 5.17 


< 1.85 






IRAS 


14059+2000 


< 0.97 


1.68 


< 1.04 


< 1.70 


< 1.11 






IRAS 


14070+0525 


< 0.53 


< 0.77 


< 1.25 


< 1.98 


< 1.37 






IRAS 


14553+1245 


< 0.89 


1.58 


< 0.89 


< 2.38 


< 1.49 






IRAS 


15327+2340 


< 1 04 


5.19 


< 24.96 


< 24.33 


< 15.22 


< 151.68 




IRAS 


16090-0139 


1.03 


3.07 


< 3.95 


1.12 


< 


2.23 






IRAS 


16255+2801 


< 5.36 


1.54 


< 12.06 


< 7.16 


< 


4.83 






IRAS 


16300+1558 


< 1.20 


0.50 


< 


1.23 


< 0.88 


< 


0.70 






IRAS 


17207-0014 


< 0.51 


6.39 


< 


6.40 


< 11.27 


< 


7.70 


12.24 


46.91 


IRAS 


18368+3549 


< 0.49 


1.18 


< 


1.07 


< 0.91 


< 


0.59 






IRAS 


18588+3517 


< 3.75 


7.96 


< 


4.91 


< 4.13 


< 


4.05 






IRAS 


20100-4156 


< 0.69 


2.86 


< 


3.58 


1.20 




1.34 






IRAS 


20286+1846 


< 0.30 


0.44 


< 


0.58 


< 0.73 


< 


0.70 






IRAS 


21077+3358 


< 0.50 


0.77 


< 


0.80 


< 0.60 


< 


0.47 






IRAS 


21272+2514 


< 0.21 


0.45 


< 


0.50 


< 0.68 


< 


0.46 






IRAS 


22055+3024 


< 0.51 


1.04 


< 


2.12 


< 2.30 


< 


1.34 






IRAS 


22116+0437 


< 1.07 


1.04 


< 


2.16 


< 2.54 


< 


2.35 






IRAS 


22491-1808 


< 0.72 


1.86 


< 


3.40 


< 7.11 


< 


3.35 


13.72 




IRAS 


23028+0725 


< 1.05 


< 1.15 


< 


1.72 


< 2.70 


< 


1.72 






IRAS 


23233+0946 


< 0.89 


3.13 


< 


1.06 


1.20 




0.66 






IRAS 


23365+3604 


< 0.44 


4.29 


< 


4.74 


< 11.96 


< 


7.97 


6.81 




IRAS 


00163-1039 


< 1.70 


30.85 


< 


7.59 


1.42 




3.35 


34.70 


73.84 


IRAS 


01572+0009 


< 0.89 


1.66 




4.76 


9.82 


< 


6.12 






IRAS 


05083+7936 


< 1.19 


18.67 


< 


1.89 


1.28 




1.33 


28.00 


71.27 


IRAS 


06538+4628 


1.17 


19.11 


< 


2.33 


< 6.97 




2.70 


37.60 


57.37 


IRAS 


08559+1053 


< 0.93 


1.35 


< 


1.34 


2.54 




0.45 






IRAS 


09437+0317 


< 0.70 


3.33 


< 


0.99 


0.47 




0.68 


9.45 


23.27 


IRAS 


10565+2448 


< 1.02 


12.42 


< 


3.07 


< 5.13 


< 


2.35 


20.89 


51.37 


IRAS 


11119+3257 


< 2.30 


< 1.75 


< 


2.40 


< 3.08 


< 


1.95 






IRAS 


13349+2438 


< 4.54 


1.68 




3.46 


7.28 


< 


5.11 






IRAS 


15001+1433 


< 0.78 


2.36 




0.67 


1.21 




0.56 






IRAS 


15206+3342 


< 1.03 


8.58 




1.36 


0.74 




1.21 






IRAS 


20460+1925 


< 0.72 


< 0.79 


< 


3.68 


1.95 


< 


1.01 






IRAS 


23007+0836 


< 3.28 


70.16 


15.51 


30.98 




9.47 


97.56 


188.36 


IRAS 


23394-0353 


< 1.10 


17.16 


< 1.92 


1.27 




2.43 


44.29 


53.48 


IRAS 


23498+2423 


< 0.55 


1.13 




1.26 


4.61 


< 


3.77 







Note. — Line fluxes are given in 10 W cm ^ . — indicates that the redshifted hne wavelength 
lay outside the range of the IRS. 
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sition at 5.51 /in. In each case, the hne number [eg, 
for H2 S(0)] indicates the rotational quantum number of 
the lower state ( J = 2 — > 0) for the quadrupole S-branch 
transition. A J — 2 results in two separate branches: or- 
tho (parallel nuclear spin, odd J) and para (anti-parallel 
nuclear spin, even J). 

We detected at least one H2 line in 49/51 OHMs and 
13/15 non-masing galaxies, with S(l) seen in all objects 
for which at least one molecular hydrogen transition is 
reported. The higher-order S(2) and S(3) lines are seen 
in roughly 2/3 of the sample, while the para ground state 
S(0) transition is detected in only ~ 15% of the sample. 
Our line detection rate is consi stent with resu l ts from the 
SINGS galaxies exami ned in iRoussel et al] (|2007D and 
the ULIRG sample of iHigdon et al.l (|2006[ ). Lower de- 
tection rates of S(0) and S(2) are likely due to a com- 
bination of the intrinsic ortho-para ratio as well as ris- 
ing continuum levels near 28 /im that can obscure weak 
line emission by S(0). IHigdon et al] (l20Q6h find that the 
S(2)/S(3) ratios are consistent with no significant differ- 
ential extinction for the two lines, which is supported 
by numerous detections in our sample of S(3) line emis- 
sion superimposed on optically deep silicate absorption 
near 9.7 /xm. We thus applied no extinction or reddening 
corrections to the line fluxes (Table [51). 

We did not detect the S(6) line in any galaxy, while 
S(4) showed a single detection in IRAS 16300-1-1558. 
Since these lines are excited by hotter gas (T 1000 K), 
they are typically weaker than the lower states which 
probe the larger reservoir of cool gas. In addition, both 
lines are only visible in the low-resolution modules at 
z ^ 0.1. This means that deblending is a significant 
issue, since both lines lie near broad PAH emission com- 
plexes. 11 OHMs and one non-masing galaxy show the 
unresolved S(7) ortho line at 5.51 /zm in the SL module. 

Measurement of the S(5) line presents a particular 
problem due to its location in a crowded section of the 
spectra. Its rest wavelength of 6.91 /im lies near the 
[Ar II] feature at 6.99 /im; in addition, both features 
are bracketed by possible hydrocarbon absorption at 6.85 
and 7.25 /tm. This not only creates difficulties in estab- 
lishing a reliable continuum, but also in deblending the 
[Ar II] and the S(5) emission (see H5.5.2I) . Emission in the 
[Ar Il]/H2 S(5) complex is seen in more than half of our 
sample, however, and so we present measurements for 
the entire feature, including blended emission from both 
lines. We caution that these fluxes should be viewed as 
upper limits for either [Ar II] or H2 S(5) emission, since 
the SL module does not have sufficient resolution to sep- 
arate the two features. 

For galaxies in which multiple H2 lines are observed, 
we fit excitation temperatur es (Tex) to the molecula r gas 
following the methods of R igopoulou et al.l ([2002') and 
IHigdon et all ([2006). Wc assume that the emission is 
optically thin (so that the lines are unsaturated), pop- 
ulations are in local thermodynamic equilibrium (LTE), 
and that the sources are unresolved in the Spitzer beam. 
The luminosity of a molecular emission line for the tran- 
sition from (J-f 2) J is then Lj = Aj x AEj x A^j+2, 
where Aj is the Einstein-A coefficient, AEj is the energy 
of the transition, and A^j+2 is the number of molecules 
in the J + 2 state. The partition function for a given 
symmetry branch is: 



X! 9 J exp[ 

Jo/p 



(2) 



where T^x is the excitation temperature and we sum only 
over a single symmetry branch (ortho or para). The sta- 
tistical weights are gj — (2 J -I- 1) x Jg, where Js = 1 for 
the para branch (even J) and Js = 3 for ortho (odd J). 

Assuming that the lines are in LTE, the ratio of level 
populations follows a Boltzmann distribution such that 
oc gj exp[—Ej/k T^x]- The inverse slope of the best- 
fit line of an excitation diagram yields T^x - Figure [5] 
shows examples of temperature fits for our data. The 
total warm H2 mass can be then calculated from T^x 
and the flux {Fj) from any transition as: 



{A^DDFjZj 



AjAEjgj exp[-Ej/kTe, 



(3) 



where (po/p is a numerical factor accounting for the ortho- 
to-para ratio (assumed to be 3:1), rn_ff2 is the mass of the 
hydrogen molecule, and is the luminosity distance. 

For cases where the S(7) line was detected, a single 
excitation temperature gives a poor fit to the full set of 
transitions. In these cases, we first fit T^x between S(3) 
and S(7), measuring hotter gas. We then subtracted this 
component from the S(0) to S(3) fluxes, and fit a second 
Ti^x to the warm gas component. This decreased the 
mean warm T^x by 20 K, with a negligible effect on 
the gas mass. We calculated the warm H2 mass using the 
flux in the S(l) transition and the hot gas mass using the 
S(3) flu x (Table El). 

Both iHigdoneLaLl (|2006l ) and IRoussel et all (|2007l ) 
suggest that the H2 emission arises from far-ultraviolet 
photons from massive stars powering photodissociation 
regions (PDRs). Detections of the H2 S(3) transition 
in nearly all objects implies that the silicate absorption 
at 9.7 /im must be partially background to the warm 
molecular gas seen in emission. Since the dust is very 
optically thick in almost all ULIRGs, this means that 
at least some molecular gas (and possibly other atomic 
transitions) actually come from superficial layers at the 
edge of the merging system. Given that the OHM is typ- 
ically formed within the central kiloparsec of the host 
galaxy, a link between the observed warm H2 gas and 
the OHM is uncertain. 

11 objects in the OHM sample and 2 non-masing 
galaxies have C O de tections publishe d in the literature 
(jSolomon et al.lll997l : lGao fc Solomonll2004al ). The beam 
width used for CO observations is several times that 
of the HR slits; since most ULIRGs are unresolved in 
the Spitzer beam, we consider the gas mass estimates 
to be comparable. The cold gas masses derived using a 
ULIRG-calibrated Mh^/Lco ratio of ~ 1.4M0/(K km 
s^^ pc^) give a warm gas mass fraction for the OHMs 
ranging from 0.04 — 0.8%, with the gas fraction of the 
non-masing galaxies lying in a similar range (0.06—0.1%). 
This is comparable to wa rm gas fractions in ULIRGs 
from IHigdon et al.l (|2006l ). implying that the mid-IR 
H2 lines probe only a small amount of the total gas mass 
in these galaxies. The bulk of the remaining portion is 
likely cold gas without sufficient energy to excite rota- 
tional transitions in the mid-IR. 
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TABLE 6 
Molecular H2 gas properties 



Object 


Ha 8(7) 


Ha S(5) 


Ha S(4) 


Ha S(3) 


Ha 8(2) 


Ha S(l) 


Ha 8(0) 




That 


^warm 


Mhot 


Xrest [Aim] 


5.51 Aim 


6.91 Aim 8.03 Aim 


9.67 fj.m 


12.28 Atm 


17.04 Aim 


28.22 fim 


[K] 


[K] 


[W Mq] 


[10^ Mq] 


IRAS 01355-1814 


< 1.52 


< 3.25 


0.22 


0.63 


1.13 


< 1.63 


262 




4.74 




IRAS 01418+1651 


< 5.29 


6.96 




1.16 


0.94 
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Note. — Line fiuxcs arc given in 10~^^ W em^^. — indicates that the redshifted line lay outside of the IRS spectral range. Fluxes for H2 S(7) 
and S(5) are measured in the SL module; all other lines are measured in the SH and LH modules. S(5) lines are tentative upper limits (indicated by 
a :) due to possible blending with [Ar 11] at 6.99 fira-. see ^5.31 
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Fig. 6. — Example H2 excitation diagrams for both non-masing (IRAS 01572+0009) and OHM galaxies (all others). Two galaxies {left) 
are fit with both warm and hot excitation temperatures; IRAS 08474+1813 fits only a warm component since the higher J lines are not 
detected. IRAS 11028+3130 shows an example of a galaxy with only a single H2 detection (for which no Tex can be determined). Dotted 
lines axe fit to the warm gas for all detections from 8(0) to S(3); dashed lines are fit to the hotter gas using detections of S(3), S(4) and 
S(7). The 8(5) line is always an upper limit due to possible blending from [Ar II] and is not used in the temperature fits. 
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5.4. PAH emission 

In addition to the atomic and simple molecular emis- 
sion lines, we also observed multiple features attributed 
to polycyclic aromatic hydrocarbons (PAHs); the broad- 
line emission com es from vibr ational modes of C-C 
and C-H bonds (jDraind l2003f) . PAH features are 
ubiquitous in the mid-IR emission of starburst galax- 
ies and U LIRGs (Lutz et al. 1998; Gcnzcl ct al. 1998; 
Sturm et al. 20 00: Peetcrs ct al. 2004; Desai et al. 2007; 
Imanishi et al.ll2007[ ). and dominate the low- resolution 
spectra of most galaxies in our sample. Multiple PAH 
features are seen for all OHMs, encompassing galaxies 
with very wide ranges in continuum shape and line emis- 
sion. We detect strong PAH transitions centered at 6.2, 
7.7, 8.6, 11.3, and 12.7 /im, several of which are also 
visible in the high-resolution spectra. Weaker emission 
features at 13.5, 14.2, 16.4, 17.1, and 17.4 /im are also 
visible in many galaxies. 

We measured the PAH emission via two methods: the 
first defines a local continuum around the PAH feature 
using a spline fit, and then integrates the total fiux af- 
ter baseline subtraction. The default continuum pivots 
are located at 5.15, 5.55, 5.95, 6.55, and 7.10 /xm for the 
6.2 /Ltm feature, and at 10.1, 10.9, 11.8, and 12.4 /xm for 
the 11.3 /im feature. These are shifted slightly for each 
object to avoid both broad absorption features (including 
water ice and hydrocarbons) and narrow atomic emis- 
sion lines. We quantify the emission from the two clean- 
est PAH features appearing in our spectra: the 6.2 and 
11.3 /im complexes (Table [T]). 

A second method uses PAHFIT (jSmith et al.ll2007D . a 
public IDL package, to fit global mid-IR spectral tem- 
plates and simultaneously measure the relative effects of 
overlapping features. The routine decomposes the low- 
resolution IRS spectra into emission from stellar contin- 
uum, dust features (including PAHs) , atomic and molec- 
ular lines, and blackbodies from thermally heated dust at 
a variety of temperatures; this is simultaneously fit with 
an extinction curve including silicate features at 9.7 and 
18 /im. Due to the large number of parameters being fit, 
however, even strong feature s may overlap suffi ciently to 
affect the accuracy of the fit ijSpoon et al.ll2002i ): in addi- 
tion, PAHFIT does not fit for several features commonly 
found in ULIRGs, such as absorption features from ice, 
hydrocarbons, and gas-phase molecules. 

PAHFIT fits the dust emission features (including 
PAHs) with a Drude profile, which has the form: 



h[X] = 



(A/A,-Ac/A)2+72 



(4) 



where Ac is the central wavelength, 7 is the fractional 
FWHM, and b is the (peak) central intensity. The Drude 
profile is typically broader than a Gaussian, with signif- 
cant amounts of power in the extended wings. Several 
dust emission features {eg, the 7.7 and 12.7 /im PAHs) 
require more than one component for a reasonable fit. 
PAHFIT returned positive detections for both the 6.2 
and 11.3 fim PAH features for nearly all galaxies; the 
fit for OHM IRAS 07572-1-0533 showed no emission at 
6.2 /im, while the fits to OHM IRAS 13218-H0552 and 
the non-masing galaxy IRAS 13349+2438 show no emis- 
sion in either dust feature. 



We compared the fiux measured in the 6.2 and 
11.3 /im PAH features from our baseline-subtracted 
spline fits to the PAHFIT values; results from PAH- 
FIT are consistently higher than those from the spline 
fit, indicating significant mixing between the PAH emis- 
sion and what was previously designated as "continuum" 
(Figure [7]). Fluxes of the 6.2 /im complex measured with 
PAHFIT are a factor of ~ 3 — 4 greater than the spline-fit 
fluxes, while the 11.3 /im feature is an average of ^ 2 — 3 
times larger. The relative strengths of the two features 
are consistent using both methods; the mean value of the 
(6.2 /im PAH/11.3 /im PAH) ratio is identical to within 
1 5% for OHMs. 

iGalliano et al.l (|2008[ ) also use both spline and multi- 
component profile fitting approaches to measure the PAH 
variations within galaxies; they show that both meth- 
ods yield the same overall trends, although each have 
their own underlying biases depending on the property 
being measured. Given the large contributions of over- 
lapping dust emission features to the 5-10 /im spectrum, 
which are not possible to separate fro m the underly- 
ing b lackbody emission using spline fits (jMarshall et al.l 
I2007f ). we consider the PAHFIT results to be the more 
robust method. Measurements of PAH features in the 
literature, ho wever, typ ically use a spline- fit method (eg , 
Brandl et al. 1120 06: Des ai et al.l [20071; iSpoon et al.l [20071; 
Zakamska et al.l 1200 8). In particular, comparisons of 



PAH data from the OHM galaxies to other s amples {eg, 
the "fork diagram" from iSpoon et al.l I2007D must use 
the same method to return physically meaningful results. 
While PAHFIT fluxes may thus better represent the ab- 
solute PAH luminosity, the spline-fit data are used when 
comparing the OHMs to objects from the literature (Ta- 
ble [7]). 

Water ice absorption at 6 /im can have significant ef- 
fects on the measurement of the 6.2 /im PAH equiva- 
lent w idth (EW); following the method of iSpoon et al.l 
(|2007f ). we correct for this by substituting the continuum 
inferred while measuring the 9.7 /im silicate strength for 
the measured 6.2 /im continuum (see tj5.5p . In total, 
24 OHMs and three non-masing galaxies with SL data 
showed absorption strong enough to affect the measured 
EW; the spline fit with the new continuum decreased the 
EW for aU objects except IRAS 11180-1-1623, for which 
the continuum levels as measured by the PAH fit and the 
silicate depth are nearly identical (within error). Since 
PAHFIT does not fit for ice absorption, we also calcu- 
late ice-corrected EW for objects showing absorption at 
6 /im by using the fiux from PAHFIT and the inferred 
continuum from silicate measurements. 

The 11.3 /im PAH is seated atop the edge of the 
deep silicate absorption at 9.7 /im; determining an 
extinction-corrected continuum level for EW measure- 
ments is thus also difficult. Spectral mapping of AGN 
galaxies with ISOCAM has shown that PAH emission 
can be spatially extended an d suppresse d near t he nu - 
cleus (jLe Floc'h et al.l 120011 ; [Dfaz-SantoseLalJ l2010[) ; 
this means that the PAH emission in AGN may be largely 
unaffected by dust absorption. Starburst galaxies, how- 
ever, can show strong PAH emission in both the 6.2 and 
11.3 /im bands in the nuclear regions where the si licate 
optical depth is at its highest ([Galliano et al.|[2008D . and 
is likely to affect continuum levels for PAH features. The 
measured 11.3 /im PAH data using the spline-fit method 
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TABLE 7 

PAH EMISSION FEATURES IN LOW-RESOLUTION SPECTRA 
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IRAS 13451+1232 


Q '^A 


y .04: 


OR 




04 
U.U4 


R 9*^ 


o. 1 u 


0.01 




01 

U.U-L 


IRAS 14059+2000 


9.16 


9.06 


0.25 




0.48 


8.93 


8.89 


0.23 




0.43 


IRAS 14070+0525 


10.37 


10.31 


1.86 


0.01 


2.41 


8.92 


9.63 


0.02 


0.01 


0.83 


IRAS 14553+1245 


y.o 1 


Q A7 


u.uo 




U.uo 


y .o ( 


9 22 


0.44 




RO 
U.uu 


IRAS 15327+2340 


10.45 


10.03 


91.96 


0.16 


8.23 


8.89 


8.67 


0.30 


0.17 


0.64 


IRAS 16090-0139 


10.17 


9.99 


1.14 


0.07 


1.45 


9.34 


9.38 


0.09 


0.07 


0.51 


IRAS 16255+2801 


9.54 


9.14 


0.94 


0.50 


1.02 


8.95 


8.58 


0.16 


0.13 


0.37 


IRAS 16300+1558 


10.31 


9.97 


1.66 


0.04 


1.09 


9.28 


9.48 


0.07 


0.04 


0.60 


IRAS 17207-0014 


10.01 


9.86 


3.58 


0.45 


2.75 


9.52 


9.24 


0.50 


0.45 


0.76 


IRAS 18368+3549 


9.82 


9.74 


2.43 




2.62 


9.50 


9.26 


0.61 




0.74 


IRAS 18588+3517 


9.84 


9.67 


1.85 


0.56 


1.65 


9.45 


9.25 


0.41 


0.23 


0.82 


IRAS 20100-4156 


9.96 


9.83 


1.52 


0.06 


1.23 


9.32 


9.38 


0.19 


0.06 


0.77 


IRAS 20286+1846 


9.44 


9.45 


1.58 


0.63 


2.26 


8.93 


8.86 


0.41 


0.20 


0.91 


IRAS 21077+3358 


9.86 


9.85 


2.88 




1.15 


9.23 


9.28 


0.21 




0.40 


IRAS 21272+2514 


9.66 


9.60 


1.97 


0.15 


1.73 


9.11 


8.97 


0.34 


0.15 


0.56 


IRAS 22055+3024 


9.22 


9.27 


0.28 




0.26 


8.88 


9.10 


0.15 




0.22 


IRAS 22116+0437 


10.27 


9.85 


3.16 




0.75 


9.34 


9.38 


0.08 




0.40 


IRAS 22491-1808 


9.47 


9.33 


1.38 


0.45 


0.98 


9.07 


8.94 


0.43 


0.45 


0.57 


IRAS 23028+0725 






















IRAS 23233+0946 


9.68 


9.59 


1.74 


0.36 


1.45 


9.30 


9.17 


0.47 


0.36 


0.76 


IRAS 23365+3604 


9.80 


9.63 


2.37 


0.27 


0.92 


9.31 


9.18 


0.35 


0.27 


0.43 


IRAS 00163-1039 


9.30 


9.16 


3.32 




0.95 


9.00 


8.83 


0.52 




0.43 


IRAS 01572+0009 


9.86 


9^47 


0.10 




0.04 


9.60 


9.55 


0.07 




0.06 


IRAS 05083+7936 


9.93 


9.95 


4.45 




2.10 


9.68 


9.57 


0.62 




0.66 


IRAS 06538+4628 


8.75 


8.77 


1.00 




0.60 


8.51 


8.50 


0.45 




0.36 


IRAS 08559+1053 


9.89 


9.86 


0.31 


0.19 


0.55 


9.63 


9.52 


0.21 


0.19 


0.29 


IRAS 09437+0317 


9.10 


9.16 


1.47 




1.99 


8.89 


8.82 


0.62 




0.77 


IRAS 10565+2448 


9.90 


9.80 


1.53 




1.43 


9.55 


9.31 


0.51 




0.51 


IRAS 11119+3257 


9.92 


9.68 


0.04 


0.05 


0.04 


9.93 


8.36 


0.06 


0.05 


0.003 


IRAS 13349+2438 












< 9.39 


< 9.10 


< 0.01 




< 0.01 


IRAS 15001+1433 


9.86 


9.86 


0.30 


0.14 


0.43 


9.54 


9.44 


0.16 


0.14 


0.20 


IRAS 15206+3342 


9.81 


9.75 


0.38 




0.41 


9.56 


9.47 


0.25 




0.27 


IRAS 20460+1925 






















IRAS 23007+0836 


9.23 


9.19 


0.33 




0.34 


8.99 


8.89 


0.19 




0.20 


IRAS 23394-0353 


9.26 


9.16 


2.24 




2.47 


8.94 


8.72 


0.53 




0.64 


IRAS 23498+2423 


9.52 


9.53 


0.06 




0.13 


9.25 


9.05 


0.04 




0.06 



Note. — The "ice" 6.2 fira PAH columns use a continuum that is corrected for water ice absorption (where present) at 6 fim. — 
indicates that PAHFIT fit no significant flux for a particular dust component. 
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IRAS 15327+2340 
I ' ' ' I ' ' 



IRAS 15327+2340 




5.6 5.8 6.0 6.2 6.4 6.6 6.8 7.0 
Wovelength (rest frome) [fim] 




5.5 6.0 6.5 7.0 

Wovelength (rest frome) [/xm] 



Fig. 7. — Differences between two methods used to measure PAH ffuxes in IRS spectra, illustrated on the spectrum of IRAS 15327+2340 
(Arp 220). (left): The observed spectrum centered on the 6.2 ^m PAH feature is shown in black. The minimum at 5.9 fj,m is the result of 
absorption by water ice. For the spline method, the PAH feature is defined as the integrated flux above a spline-interpolated continuum 
{green dotted line) between 5.95 and 6.5 fim. The dot-dashed blue lines represent individual Drude fit components from PAHFIT, while 
the red line is the continuum emission (blackbody dust -|- starlight). The dashed grey line shows PAHFIT's global fit to the data, (right): 
PAHFIT (blue) and spline-fit (green) results with the continuum subtracted; the feature as measured by PAHFIT has significantly more 
flux than with the spline continuum, typical of nearly all galaxies in our sample. 
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are therefore hkely to underestimate the himinosities. 

5.5. Absorption features 
5.5.1. Silicates 

The LR spectra show near-ubiquitous absorption 
from amorphous sihcate dust, with a strong feature 
caused by a Si-0 stretching mode near 9.7 /im and 
a weaker fea ture caused by an Si-O-S i bending mode 
near 18 /xm ()Knacke fc ThomsonI HItI . The presence 
of dust is unsurprising, as the characteristic extreme IR 
luminosities of ULIRGs are caused by large a mounts 
of hea ted dust being thermally re-radiated. Ha o et ahl 
()2007t ) found that ULIRGs nearly uniformly show ab- 
sorption in the two silicate features, in contrast to QSOs 
and some Seyfert galaxies which typically show the fea- 
ture in emission ( Siebenmorgcn et al. 2005; JI3ao et ahl 
120051: iSturm et aLi r2005: Schweitzer et al. 200i). 



We measure the strength of the silicate ab sorption at 
both 9.7 and 18 ^m using the method of iSpoon et al.l 
(|200l : 



fx 



fci 



(5) 



where f\ is the measured flux and fcont the interpolated 
continuum at the feature extremum. More negative val- 
ues of Ssii represent deeper absorption. The expected 
continuum is calculated using a combination of spline and 
power-law fits, depending on the str ength of the PAH a nd 
water ice features in the spectrum (jSpoon et al.ll2007[ ). 

All OHMs and ^ 95% of the non-masing galaxies 
showed absorption at both 9.7 and 18 /im (Table H]); 
the average depth for the OHMs is S'9.7 = —1.8 ± 0.8, 
while the average depth of the non-masing galaxies is 
5*9 7 = —0.6 ±0.4. The deepest absorption is in the 
OHM IRAS 04454-0838 (S'9.7 = -3.7), while only one 
object (the non-masing galaxy IRAS 13349-1-2438) shows 
emission in both amorphous silicate features. 

In addition to the amorphous silicate, we detect weaker 
features from crystalline silicate absorption in 19 OHMs, 
including bands at 11, 16, 19, 23, and 28 / im (see Fig- 
ure El f or an example). We use the method of'S poon et al.l 
()2006l ) to subtract off both the dust continuum and 
amorphous component to measure the residual optical 
depth at 16 and 23 fim, typically the strongest crys- 
talline features (Table [S]). The deepest Siq occurs for 
IRAS 10039-3338, at -0.4; however, the S/N ratio 
means we are only sensitive to an absorption limit of 
5*16 — —0.1 and Sis — —0.05. Only one detection of 
crystalline silicates is made in a non-masing galaxy, in 
IRAS 11119+3257. Since all detections of crystalline sil- 
icates in OHMs have 5'9.7 < —1.4, the lower total dust 
column in non-masing galaxies is a likely contributor to 
the detection rate. 

5.5.2. Aliphatic hydrocarbons 

Absorption bands arising from hydrogenated amor- 
phous carbon grains (HACs) can also be sig nificant con- 
tributors to diff use dust in the gala c tic ISM (jChiar et al.l 
[200m. ISO (ISpoon et al.l [2OOIL [200I a nd Sv itzer 
partois et all l2007t iDartois fc Mufioz-Carol [20071) ob- 
servations have identified HAG absorption features 
due to bending modes at 6.85 /Ltm (GH2/GH3) and 



7.25 (GH3) in ULIRGs. These aliphatic features 
represent a counterpart to the aromatic hydrocarbons 
responsible for PAH emission and are an abundant com- 
ponent of the ISM in luminous galaxies. 

We detect absorption from the 6.85 /im HAG transi- 
tion in 27/51 galaxies in the OHM sample (Figure H]), 
with zero detections in the non-masing sample. The 
accopmanying 7.25 /tm feature is detected in eight of 
the galaxies in which the 6.85 /tm feature is seen. The 
strength of the HAG is measured using a spline fit with 
pivots at 5.2, 5.6, 7.8, 14.0, and 26.0 /im to determi ne the 
local continuum in the 6-8 /im region (S poon et af]l2007[ ) 
and integrate the total fiux within the absorption feature 
(TablelH). The average optical depth of the 6.85 /tm fea- 
ture is T6.85 = 0.23 ± 0.11, and the average depth of 
the 7.25 /tm feature is r7.25 = 0.20 ± 0.06. Many of the 
galaxies with no detection of HAGs, however, have limits 
from noise that are consistent with the absorption depths 
measured in brighter galaxies. 

5.5.3. Ices 

Absorption from ices in a variety of molecular species 
(including H2O, GO, GO2, and GH4) has been detected 
in spectroscopy of IR-br ight galaxies (Sp oon et al.ll200dl 
120011: ISturm et al]|2000[ ). The band from water ice ab- 
sorption stretching from 6-8 /tm is prominent (Figure [5]) 
and was detected in ^ 10% of a sample of b right galaxies 
using ISO ISpoon et al. 20021) and the IRS (lArmus et al.l 
12004; Spoon et al. 2005: lArmus et al.ll2007l ). 

We detected water ice absorption at 6 /tm in 24 OHMs 
and three non-masing galaxies. We use the spline contin- 
uum from fitting the 9.7 /tm silicate feature as the local 
5.5-7 /tm continuum in order to obtain an optical depth 
spectrum for the 6 /im absorption complex. The result- 
ing water ice optical depths are tabulated in Table El We 
note that contamination by 6.2 /tm PAH emission and 
absorption by other species than water ice (|Spoon et al.l 
l2005h may add confusion in properly measuring optical 
depths. 

5.5.4. C2H2, HCN, and CO2 

Previous mid-IR surveys have also identified bands of 
molecular ga s abso r ption in ULIRGs (ISpoon et al.l 120061 : 
lArmus et all 120071: ILahuis et all I2007D . including the 
vibration-rotation bands of acetylene (G2H2; 13.7 /tm), 
hydrogen cyanide (H GN; 14.02 /xm ) , and carbon diox- 
ide (GO2; 15.0 /im). ILahuis et all ( |200'<t ) reported the 
detection of both G2H2 and HGN in fifteen (U)LIRG nu- 
clei, with detections of GO2 in four objects. Eight of 
the objects in the Lahuis sample are OHMs in our sam- 
ple; we confirm detections of G2H2 in all eight galax- 
ies, in addition to the OHMs IRAS 10039-3338 and 
IRAS 12018-fl941 (Figure El)- HGN is detected in only 
4/10 archival galaxies, meaning that we cannot confirm 
the HGN detection of four galaxies; since the optical 
depth of HGN is typically much weaker than that of 
G2H2, however, it is possible that our lower detection rate 
is a result of improved S/N in their reduction process. We 
also confirm the detection of GO2 in IRAS 15327-f 2340 
( Arp 220) . No galaxies in the non-masing control sample 
showed absorption in any molecular band, nor did any of 
the OHMs observed in our dedicated program. 

Since these gas-phase absorption features are actually 
a blend of multiple absorption lines, the peak optical 
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TABLE 8 
Solid-phase absorption features 



Object 


6.0 iim H2O ice 


6.85 


^tm HAG 7.25 


fim HAG 


9.7 Aim 


18 fim 


16 /im 


23 fim 




r 


r 


flux T 


flux 






Qresid 
sil 


Qresid 
sil 


IRAS 01355-1814 










-2.4 


-0.9 


-0.3 


-0.2 


IRAS 01418+1651 




0.49 


-8.3 




-1.3 


-0.4 






IRAS 01562+2528 










-0.7 


-0.3 






IRAS 02524+2046 




0.20 


-0.4 




-0.9 


-0.3 






IRAS 03521+0028 


0.42 


0.13 


-0.4 




-1.4 


-0.2 


-0.1 


-0.1 


IRAS 04121+0223 


1.61 








-1.0 


-0.2 






IRAS 04454-4838 


0.42 


0.35 


-9.0 0.10 


-1.5 


-3.7 


-1.0 


-0.3 


-0.2 


IRAS 06487+2208 




0.23 


-4.6 




-1.2 


-0.3 






IRAS 07163+0817 










-1.2 


-0.1 






IRAS 07572+0533 










-0.6 


-0.3 






IRAS 08201+2801 


1.06 


0.37 


-4.0 0.22 


-1.5 


-2.2 


-0.6 


-0.2 


-0.1 


IRAS 08449+2332 










-1.2 


-0.5 


-0.1 


-0.1 


IRAS 08474+1813 




0.36 


-0.2 




-1.9 


-1.2 






IRAS 09039+0503 


0.98 


0.15 


-0.5 




-2.0 


-0.6 






IRAS 09539+0857 




0.24 


-2.5 




-3.1 


-1.2 


-0.4 


-0.2 


IRAS 10035+2740 










-1.5 


-0.8 






IRAS 10039-3338 


0.23 


0.23 


-166.3 




-3.1 


-1.0 


-0.4 


-0.3 


IRAS 10173+0828 










-1.9 


-0.8 


-0.3 


-0.2 


IRAS 10339+1548 










-1.1 


-0.05 






IRAS 10378+1109 


0.72 


0.18 


-0.6 




-2.0 


-0.3 






IRAS 10485-1447 




0.25 


-0.4 




-2.9 


-0.9 






IRAS 11028+3130 










-2.6 


-1.0 






IRAS 11180+1623 


0.54 








-1.7 


-0.5 






IRAS 11524+1058 




0.27 


-1.5 0.2: 


-0.3: 


-1.5 


-0.8 






IRAS 12018+1941 




0.16 


-2.6 




-1.4 


-0.4 






IRAS 12032+1707 


0.71 


0.55 


-6.6 0.30 


-3.9 


-2.7 


-0.8 






IRAS 12112+0305 


0.59 


0.41 


-3.5 




-1.8 


-0.3 


-0.2 


-0.1 


IRAS 12540+5708 










-0.7 


-0.2 






IRAS 13218+0552 










-0.5 


-0.4 






IRAS 13428+5608 


0.50 


0.40 


-28.5 




-2.0 


-0.5 






IRAS 13451+1232 










-0.5 


-0.1 






IRAS 14059+2000 










-0.8 


-0.1 






IRAS 14070+0525 


0.90 


0.24 


-1.8 0.15 


-1.8 


-2.7 


-0.9 






IRAS 14553+1245 










-1.3 


-0.5 






IRAS 15327+2340 


0.68 


0.35 


-50.1 




-3.1 


-0.4 


-0.2 


-0.1 


IRAS 16090-0139 


0.56 


0.45 


-10.5 0.24 


-5.6 


-2.4 


-0.6 






IRAS 16255+2801 


0.54 








-2.2 


-0.6 






IRAS 16300+1558 


0.61 


0.38 


-2.2 0.21 


-1.1 


-2.7 


-0.7 


-0.3 


-0.1 


IRAS 17207—0014 


0.31 


0.23 


-17.4 




-1.9 


-0.6 


-0.2 


-0.1 


IRAS 18368+3549 










-1.8 


-0.2 


-0.3 


-0.2 


IRAS 18588+3517 


0.72 








-2.2 


-0.6 


-0.3 


- 


IRAS 20100-4156 


1.45 


0.23 


-2.9 0.23 


-4.5 


-2.4 


-0.7 


-0.2 


-0.1 


IRAS 20286+1846 


1.08 








-1.6 


-0.6 






IRAS 21077+3358 










-1.9 


-0.7 






IRAS 21272+2514 


1.66 


0.21 


-0.7 




-2.8 


-0.7 






IRAS 22055+3024 










-1.3 


-0.3 






IRAS 22116+0437 




0.04; 


-0.1: 




-2.6 


-0.9 


-0.2 


-0.2 


IRAS 22491-1808 


0.43 


0.19 


-0.2 




-1.5 


-0.5 


-0.2 


- 


IRAS 23028+0725 












-0.6 






IRAS 23233+0946 


0.44 








-1.9 


-0.4 






IRAS 23365+3604 


0.66 








-2.0 


-0.5 






IRAS 00163-1039 










-0.5 


-0.1 






IRAS 01572+0009 










-0.2 


-0.2 






IRAS 05083+7936 










-1.1 


-0.3 


















n 9 






IRAS 08559+1053 


0.18 








-0.6 


-0.2 






IRAS 09437+0317 










-1.1 


-0.3 






IRAS 10565+2448 










-1.2 


-0.3 






IRAS 11119+3257 


0.19 








-0.7 


-0.3 


-0.2 




IRAS 13349+2438 










0.1 


0.07 






IRAS 15001+1433 


0.30 








-0.9 


-0.4 






IRAS 15206+3342 










-0.4 


-0.2 






IRAS 20460+1925 












-0.4 






IRAS 23007+0836 










-0.3 


-0.1 






IRAS 23394-0353 










-0.7 


-0.3 






IRAS 23498+2423 










-0.6 


-0.3 






Note. — Silicate s 


trcngth is defined in 


Eqn. [5] the 9.7 and 18 


/^ni features 


are depth 


LS for aniorpiious 


silicates, 



while the 16 and 23 fim crystaUiiie features arc residual depths measured after the 18 fim feature was subtracted. 
Fluxes are given in cm ; objects marked with a : represent uncertain detections. 
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30.0 



33.5 34.0 34.5 35.0 14.0 

Wavelength [fj,m] 



16.0 



18.0 



20.0 



Fig. 8.— Examples of mid-IR absorption features for the OHM IRAS 15327+2340 (Arp 220). Top left: The 6.0 fim H2O ice and 
6.85 fj,m HAG absorption features. Top right: Gas-phase G2H2 13.7 /im and HCN 14.0 fim absorption. Bottom left: Gas-phase OH 
34.6 fim absorption. Bottom right: Residual optical depth of the 16 p.m crystalline silicate feature. Measurements for all labeled features 
are in Tables [S] and [9] 
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TABLE 9 

Gas-phase absorption features 



Object 


13.7 /im 


C2H2 


14.02 ura 


HCN 


15.0 Atm CO2 


fnorm 


flux 


fnorm 


flux 


fnorm flux 


IHAb 10039—3338 


0.93 


—0.9 


0.86 


-4.0 




IRAS 12018+1941 


0.95 


-0.7 


0.94 


-0.8 




IRAS 12540+5708 


0.95 


-9.5 








IRAS 13218+0552 


0.82 


-1.1 








IRAS 13428+5608 


0.93 


-0.8 








IRAS 14070+0525 


0.79 


-0.4 








IRAS 15327+2340 


0.84 


-5.8 


0.90 


-5.6 


0.94 -2.0 


IRAS 16090-0139 


0.86 


-1.1 








IRAS 17207-0014 


0.93 


-0.7 








IRAS 20100-4156 


0.82 


-1.4 


0.84 


-1.1 




IRAS 22491-1808 


0.92 


-0.7 








Note. — fnorm gives the peak depth of absorption features plotted in nor- 



1 1 



10 
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T3 



L C2H2 


' ' 1 
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CO2 J 




. . 1 
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13.0 13.5 14.0 14.5 15.0 
Wavelength (rest frame) [/^m] 



15.5 



malizcd fiux units. Fluxes arc measured in 10 



W em 



depth measured is a function of the velocity resolution 
of the spectrograph. We therefore report the integrated 
fl-ux and the peak depth in normalized flux units {fnorm, 
where the spectrum has been divided by the adopted 
continuum; Spoon ct al. 2004) in Table [9l 

iLahuis et al.l (|2007D model abundances for ULIRGs 
with detections of C2H2, HCN, and CO2, and suggest 
that they are associated with a phase of deeply embed- 
ded star formation, excluding the possibility of the fea- 
tures arising from an X-ray dom inated region (XDR) 
powered by AGN. iDarlingI (J2007) has also shown that 
OHMs have the highest mean molecular gas densities 
among starburst galaxies (traced by the J = 1 — > 
rotational HCN transition) and also possess high dense 
molecular gas fractions, comprising a dis tinct popula- 
tion in the IR-CO relation. The results of ILahuis et al.l 
(jlOO?) show that 9/15 ULIRGs with absorption in both 
C2H2 and HCN are known OHMs in our Spitzer sam- 
ple; this dense gas fraction (~ 50%) is nearly identical 
to the observed OHM fraction in starbursts with dense 
(ULIR Gs with Lhcn I Lnn > 0.07) f ractions of m olecu- 
lar ga s (|Gao fc Solomonll2004bl : lDarlin g 2007; Baau'eTaD 
[2008^ . 

Given that the S/N ratio for the OHM and non-masing 
galaxies are of comparable magnitude, the lack of detec- 
tion of any gas-phase species in the non-masing galax- 
ies is a striking difference compared to the OHMs. Fig- 
ure ini shows the median stack of both samples near the 
regions of gas-phase absorption; while the C2H2 feature 
at 13.7 /^m can be clearly seen in the median OHM spec- 
trum, neither the HCN nor the CO2 transition is promi- 
nent. Since the data have been median stacked (as op- 
posed to a mean, which can be dominated by a few deep 
absorbers), this suggests low-levels of C2H2 present in a 
significant fraction of the OHM host galaxies. No molec- 
ular absorption appears in the medianed spectrum for 
the non-masing galaxies. 

While a connection between OHMs and dense molecu- 
lar gas is known to exist, the lack of detected molecular 
absorption in the mid-IR for the majority of OHMs is not 
entirely unexpected. OHMs occur in merging galaxies 
where different populations of gas may be kinematically 
and thermally distinct, yet are observed a s a single un- 
resolv ed region within the Spitzer beam. ILahuis et al.l 
(|2007f) suggest that high abundances of warm, dense 
gas are associated with deeply-embedded star formation, 



Fig. 9. — Medianed HR spectra for both OHMs {black) and non- 
masing {grey) galaxies, sampled at intervals of 0.01 /xm and nor- 
malized in flux at 15 lira. PAH emission is visible in bands cen- 
tered at 13.6 and 14.2 fira. The dotted lines mark locations of 
gas-phase absorption in C2H2 (13.7 lira), HCN (14.02 /im), and 
CO2 (15.0 /^m). The spectra are vertically offset to highlight the 
differences between the samples. 

where HII regions are prevented from expanding by large 
pressure gradients and extend the lifetime of the star for- 

mation proce ss. 

iBaan et al] ([2008) interpret dense gas abundances as 
excluding very hard radiation fields (such as those found 
in XDRs) that dissociate the molecules; they suggest that 
the molecular emission arises from PDRs surrounding 
HII regions. Although few OHMs show absorption from 
dense molecular gas, the OHM sample also has few iden- 
tified AGN or XDRs (only 4/51 OHMs show [Ne V] at 
14 /im). The connection between dense molecular gas 
and the presence of an AGN is thus unclear based on 
this data alone. 



2n 



5.5.5. Gas-phase OH 

For three OHMs, we report detection of the 
J=5/2 2^3/2 J=3/2 OH absorption doublet 
near 34.616 /zm: III Zw 35 (IRAS 01418-^1651), Mrk 273 
(IRAS 13428+5608), and Arp 220 (IRAS 15327-^2340). 
This feature is generally difficult to detect since it lies 
near the noisy, far-red edge of the LH module; for all ob- 
jects with z > 0.08, it is redshifted out of the IRS range. 
O H absorption at 34.6 u m in Arp 220 was first reported 
bv iSkinner et all (119971) u sing I SO and confirmed with 
the IRS bv IParrah et al.] (pOOl . who incorrectly iden- 
tified it as the 0H~ ion. All OH absorption features 
are well fit with a single Gaussian, since the separation 
between the doublet features (AA ~ 0.02 fj,m) is com- 
parable to the resolution element in the LH module. No 
detection of OH absorption was made for any of the non- 
masing galaxies. 

Assuming the OH transitions are optically thin, we can 
use the EW to derive a column density for the OH ground 
state, which is likely to be a good proxy f or the total col- 
umn at typical molecular cloud densities (jBradford et al.l 
IT999h : 



EW Sttc gi 



(6) 



Aui A4 g„" 

OH column densities for all galaxies are quite similar, 
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TABLE 10 

Properties of OH gas-phase absorption 



OHM foH Tp^ak EW NoH labs lOHM (ppump 

[10-21 W cm-2] [10-3 fim] [cm-2] [ph/s] [ph/s] % 

IRAS 01418+1651 (III Zw 35) -6.1 0.15 5.5 1.1 x 10^^ 1.7 x 10^4 1.8 x 10^^ 10 

IRAS 13428+5608 (Mrk 273) -28.5 0.14 10.2 2.1 x lO^^^ 1.7 x 10^5 1.4 x 10^^ 0.8 

IRAS 15327+2340 (Arp 220) -172 0.21 15.0 3.0 x 10^'^ 2.3 x 10^^ ^ g ^ lO''^ 0.7 



Note. — The pumping efficiency {4>pump — Ioh m Mats x 100) assumes all pumping comes from the 34.6 /xm transition. 



lying between (1 - 3) x lO" cm"^ (Table [TU]). The 
measured Nqh from the IRS data for Arp 220 also 
agree s within a factor of 2 of the column measured with 
ISO (jSkinner et al.|[T997[) . Limits for galaxies in which 
the 34.6 fim OH feature is not detected are of order 
NoH < 1 X 10" cm-2. 

We compare the Nqh derived from the rotational 
34.6 /.tm transitions to the OH column density measured 
in galaxies who show the hyperfine 1667 MHz feature in 
absorption. The majority of such galaxies are ULIRGs 
of comparable luminosity to the galaxies in our non- 
masing sample. Measurements from ten OH absorbers 
(jBaan et al.l(l99l iDarlingI [20071 ) give Nqh = rex(l-8 ± 
1.9) X 10^^ cm-2, ^j^gj-g j"^^ jg lY^Q Qjj excitation temper- 
ature in K. If the dust and gas are well-mixed, then the 
temperature of the dust (~ 50 — 100 K) can be used as a 
proxy for T^x- This gives OH column densities for both 
OHMs (34.6 //m) and non-masing galaxies (1667 MHz) 
with comparable values of Nqh — 10^^ cm^^ jf go, 
then this addresses one of the crucial differences between 
OHMs and non-masing ULIRGs - namely, that differ- 
ences in the abundance of masing molecules are not a 
key factor for triggering an OHM. 

The amount of OH available in the galax y can also test 
models of the OHM pumping mechanism. iSkinner et alJ 
(|1997f ) computed the photon flux (7abs = L^b^ /hvon) 
absorbed in the 34.6 transition from the OHM 
Arp 220. They found that jabs is roughly 1% of the pho- 
ton flux in the OHM {'Johm ~ Lqhm /hvoHAi)- If the 
18-cm and mid-IR pumping photons lie along the same 
line of sight, this means that pumping photons from the 
34.6 /im transition alone can power the OHM (given an 
efficiency of - ^^ 1% or higher). While radiative transfer 
models from iLockett fc Elitzuii ()2008[ ) suggest that the 
53 /im OH transition likely contributes more pumping 
photons than the 34.6 fim line, the energetics are consis- 
tent with the basic acccepted mid-IR pumping model. 

For the galaxies in our sample with OH absorption, 
Arp 220 and Mrk 273 would require pumping efficien- 
cies on the order of 1% to power the OHMs from the 
34.6 fim transition alone. HI Zw 35 shows the weakest 
OH absorption among our detections and would require 
an efficiency of (f>pump — 10%. 

6. CONCLUSIONS 

We present mid-infrared spectra and photometry for 
51 OH megamasers taken with the IRS on Spitzer, along 
with 15 galaxies confirmed to have no megamaser emis- 
sion above Lqh — IO^-^Xq. All objects in both samples 
have full coverage in both the low- and high-resolution 
IRS modules. We measure both emission (PAH, H2, 
and fine-structure atomic transitions) and absorption 



(silicates, hydrogenated amorphous carbon grains, and 
molecular bands) features, with full spectra, line fluxes, 
equivalent widths, and absorption depths presented for 
each object. 

The majority of the galaxies closely resemble standard 
mid-IR ULIRG templates, with the low-resolution emis- 
sion dominated by moderate-to-deep amorphous silicate 
absorption at 9.7 and 18 /im and PAH features at 6.2, 
7.7, 8.6, 11.3 and 12.7 /im. The OHMs (on average) 
show deeper silicate absorption and steeper continuum 
slopes than the non-masing galaxies. Crystalline sili- 
cate absorption is detected in roughly a third of OHMs, 
but in only 1 out of 15 non-masing galaxies. OHMs are 
also the only galaxies in our sample to show absorption 
from hydrogenated amorphous carbon (HAG), gas-phase 
HGN, G2H2, and GO2; however, higher average noise 
in the non-masing spectra mean that features at simi- 
lar absorption depths could be obscured due to lack of 
sensitivity. High-resolution spectra show emission from 
[Ne II] and [Ne III] in almost all galaxies, with emission 
from [S III], [S IV], and [O IV] commonly detected. The 
high-ionization [Ne V] line (a clear tracer of AGN) is de- 
tected in < 10% of OHMs and in 53% of the non-masing 
galaxies. Almost all galaxies in both samples also show 
emission in multiple H2 rotational transitions. 

We also measure the 34.6 /im OH transition in three 
OHMs. OH column densities derived from the mid-IR 
OH transition are of the same order of magnitude as the 
column densities derived from the 1667 MHz OH tran- 
sition for ULIRGs in the literature. We interpret this 
as evidence that the OH abundances in both OHMs and 
non-masing galaxies are similar, and are not a limiting 
factor for megamaser emission. 

A companion paper (Willett et al. 2011; Paper II) 
presents a full analysis of the mid-IR data, with compar- 
isons between the two samples and connections to the 
OH megamaser properties. 

This work is based on observations made with the 
Spitzer Space Telescope, which is operated by the Jet 
Propulsion Laboratory, Galifornia Institute of Technol- 
ogy under a contract with NASA. Support for this 
work was provided by NASA through grant 30407 is- 
sued by JPL/Galtech. We have made extensive use of the 
NASA/IPAG Extragalactic Database (NED) which is op- 
erated by JPL and Galtech under contract with NASA. 
Many thanks are due to J.-D. Smith for help with PAH- 
FIT, N. Halverson for comments on computing flux lim- 
its, D. Farrah for useful discussions, and to the Spitzer 
Science Genter for hosting KWW and JD in April 2007. 
VG acknowledges partial support from the EU ToK grant 
39965 and FP7-REGP0T 206469. 
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Fig. 10. — The 34.6 OH absorption feature in III Zw 35, Mrk 273, and Arp 220. All spectra are been normalized in flux near 34.5 ^m. 



APPENDIX 

HIGH-RESOLUTION DATA WITH BACKGROUND SKY SUBTRACTION 

As discussed in 331 the reduction process for the overah sample is shghtly different for some archival galaxies that did 
not have separate IRS sky backgrounds in the high-resolution modules. Since much of our subsequent analysis (Paper 
II) depends on statistical comparisons between the two samples, we chose to minimize possible systematic errors and 
reduced all galaxies in a uniform manner without HR sky subtraction. These data are, however, likely to be a more 
reliable indicator of the absolute flux levels due to subtraction of the zodiacal background; therefore, we also present 
atomic and molecular line fluxes for these galaxies. 
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TABLE 11 

Hl-RES LINE FLUXES FOR COMMON ATOMIC EMISSION LINES WITH HR SKY 
SUBTRACTION 



Ohiprt 

\J UJ Cly U 


[O iVJ 


[Nc ii] 


[Nc III] 


fS tttI 
[O iiij 




10 511 iim 




1 ^ / /TTl 
A-^ttJ^rJ Ll/Lll 


1 8 71 3 iim 


IRAS 01562+2528 


— 


— 


— 


0.65 


IRAS 02524+2046 


— 


— 


— 


0.47 


IRAS 04121+0223 


- 


- 


- 


0.67 


IRAS 04454-4838 


0.42 


1.95 


0.43 


_ 


IRAS 06487+2208 


_ 


_ 


_ 


5.01 


IRAS 07163+0817 


— 


— 


— 


1.34 


IRAS 08201+2801 




_ 


_ 


0.41 


IRAS 08449+2332 


_ 


_ 


_ 


1.60 


IRAS 08474+1813 


_ 


_ 


_ 


0.14 


IRAS 10035+2740 


_ 


_ 


_ 


0.35 


IRAS 10039-3338 


0.98 


17.22 


4.20 


8.08 


IRAS 10339+1548 


- 


- 


- 


0.97 


IRAS 11180+1623 


- 


- 


- 


0.40 


IRAS 11524+1058 


- 


- 


- 


0.27 


IRAS 12540+5708 


— 


19.47 


— 


— 


IRAS 14059+2000 








0.59 


IRAS 14553+1245 








1.58 


IRAS 15327+2340 




59.39 


6.73 


7.54 


IRAS 16255+2801 








1.54 


IRAS 18368+3549 








1.01 


IRAS 18588+3517 








2.89 


IRAS 20286+1846 








0.48 


IRAS 21077+3358 








0.74 


IRAS 21272+2514 




2.22 


0.33 


0.40 


IRAS 22055+3024 








1.16 


IRAS 22116+0437 








0.90 


IRAS 00163-1039 


■>.l, 1 


87.43 


14.30 


32.10 


IRAS 05083+7936 




49.95 


7.63 


19.60 


IRAS 06538+4628 


0.84 


47.34 


5.90 


20.49 


IRAS 09437+0317 




8.72 


1.22 


3.58 


IRAS 23394-0353 




46.45 


7.60 


18.50 



Note. — Fluxes are in cm . No data is given for the 26 OHMs in 

our program for lines with Xrest ^ 16 jum since we have no SH sky subtraction 
available for these objects. 



TABLE 12 

Hl-RES LINE FLUXES FOR RARER ATOMIC EMISSION LINES WITH HR SKY SUBTRACTION 



Object 


H I 7-6 


[Ne v] 


[01 n] 


[Fe II] 


[Ne v] 


[O IV] 


[Fe n] 


[S in] 


[Si II] 


Xreat [MIIi] 


12.368 


14.322 


14.369 


17.936 


24.318 


25.890 


25.988 


33.481 


34.815 


IRAS 04454-4838 












1.12 




0.90 




IRAS 06487+2208 












1.12 








IRAS 07163+0817 








0.73 












IRAS 10339+1548 










0.84 


2.55 








IRAS 11524+1058 








0.35 




0.59 








IRAS 14553+1245 








0.41 












IRAS 16255+2801 












0.46 








IRAS 18368+3549 








0.51 












IRAS 21272+2514 


0.20 


0.10 
















IRAS 00163-1039 


0.64 










1.56 


2.46 


32.18 


74.28 


IRAS 05083+7936 




0.58 


0.43 






1.53 


1.97 


29.34 


43.28 


IRAS 06538+4628 




1.08 




1.06 






2.15 


37.21 


48.72 


IRAS 09437+0317 












0.47 


0.68 


9.60 


22.74 


IRAS 23394-0353 


0.59 










1.49 


2.53 


42.37 


53.56 



Note. — Fluxes are in cm . No data is given for the 26 OHMs in our program for lines with 

Xrest ^ 16 jum since we have no SH sky subtraction available for these objects. 
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TABLE 13 

Hl-RES LINE FLUXES AND UPPER LIMITS FOR, H2 TRANSITIONS 
WITH HR SKY SUBTRACTION 



Object 


H2 S(3) 


H2 S(2) 


H2 S(l) 


H2 S(0) 


\reat [/Um] 


9.67 


12.28 


17.04 


28.22 


IRAS 01562+2528 






0.78 




IRAS 02524+2046 






0.72 




IRAS 04454-4838 


1.05 


1.21 


3.03 


0.76 


IRAS 06487+2208 






2.03 




IRAS 08201+2801 






0.51 




IRAS 08449+2332 






1.40 




IRAS 08474+1813 






0.27 




IRAS 10035+2740 






1.16 




IRAS 10039-3338 


3.39 


1.79 


3.85 




IRAS 10339+1548 






0.46 




IRAS 11180+1623 






0.89 




IRAS 11524+1058 






0.81 




IRAS 12540+5708 


2.42 


4.22 






IRAS 14059+2000 






2.55 




IRAS 15327+2340 




7.36 


15.42 


14.55 


IRAS 16255+2801 






0.55 




IRAS 21077+3358 






1.17 




IRAS 21272+2514 


0.26 


0.35 


0.68 




IRAS 22116+0437 






1.35 




IRAS 23028+0725 






1.11 




IRAS 00163-1039 


2.18 


2.82 


6.01 




IRAS 05083+7936 


2.34 


2.64 


5.11 




IRAS 06538+4628 




3.33 


8.70 


2.99 


IRAS 09437+0317 






2.67 


1.78 


IRAS 23394-0353 


3.40 


2.25 


4.99 


1.38 



Note. — Fluxes are in 10~^^ W cm~^. No data is given for the 
26 OHMs in our program for the S(2) or S(3) Unes since we have 
no SH sky subtraction available for these objects. 
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